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Einleitung 



In den letzten Jahren hatte sich herausgestellt, dass zwischen den Bose-Einstein-Kondensaten 
(BEC) und Objekten aus der Allgemeinen Relativitatstheorie (ART) eine Analogie besteht 
|17j . Die Bose-Einstein-Kondensation konnte erstmals von Anderson et al. (1995) ^ experimental! 
iiberpriift werden. Seitdem sind in Experiment und Theorie groBe Fortschritte zu verbuchen. 
Die Voraussetzungen, Kondensate in vielen verschiedenen Konfigurationen herzustellen, wurden 
mdglich. Anhand dreier solcher Kondensate - zigarrenformig, ringformig und frei expandierend 
- wird in dieser Arbeit die Verkniipfung zur Gravitationsphysik gezeigt. Es konnen zwei 
Eigenschaften im BEC, die aus der Physik in gravitativen Feldern bekannt sind, im Kondensat 
entdeckt werden. Die Bewegungsgleichung einer von auBen angeregten Storung kann im 
hydrodynamischen Limes mit einer efFektiven Metrik formuliert werden. Abhangig vom 
gewahlten Kondensat kann die Metrik einem Objekt aus der Allgemeinen Relativitatstheorie 
zugeordnet werden, z.B. einem Scliwarzen Loch und dem de-Sitter-Universum. AuBerdem 
kann gezeigt werden, dass in Anwesenheit dynamischer Instabilitaten Moden im Kondensat 
angeregt werden. Es handelt sich dabei um Quasiteilchen, den Phononen, die paarweise mit 
positiver und negativer Energie erzeugt werden. 

In den ersten beiden Kapiteln werden die Grundlagen zur Bose-Einstein-Kondensation 
und der Allgemeinen Relativitatstheorie behandelt. Der Hauptteil beginnt mit Kapitel (jlllj) . 
der Verifizierung des Zusammenhangs zwischen BEC und Gravitation. Im Anschluss daran 
werden drei spezielle Kondensate analysiert. In Kapitel ()IV|) wird ftir das zigarrenformige 
und das ringformige Kondensat die Metrik berechnet und mit der eines Schwarzen Lochs 
verglichen. In beiden Systemen wird gezeigt, dass Instabilitaten auftauchen und wie diese sich 
mit der Zeit entwickeln. In Kapitel (|V|) ist das System ein frei expandierendes Kondensat. Es 
wird eine Metrik gesucht, die das de-Sitter-Universum beschreibt. Am Ende der Kapitel ()IVI) 



und fyj) werden die Ergebnisse jeweils zusammengefasst. Den Schlufipunkt setzt eine kurze 
Diskussion aller erzielten Resultate und ein Ausblick darauf, was aktuell in der Forschung 
von Inter esse ist. 
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Teil I 

Einfuhrung in die 
Bose-Einstein-Kondensation 
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Kapitel 1 



Das ideale Bose-Gas 



Als ideal bezeichnet man ein Gas, bei welchem die Atome nicht miteinander wechselwirken. 
Nehmen wir zusatzlich an, die Atome sind in einem auBeren Potential, so konnen sie nur 
bestimmte Energiewerte annehmen (Abb. [TTT|) . 

Von Interesse ist nun die Frage, wie die verschiedenen Energieniveaus bei gegebener 



Energie 




Abbildung 1.1: Beim Ubergang von der klassischen Beschreibung (a) eines Gases in den 
quantenmcchanischen Formalismus treten zwei Sorten von Teilchen auf. Wahrend die Bosonen^ (b) 
jeden Zustand beliebig oft besetzen konnen. ist cs zwei Fermioncn^ (c) nicht erlaubt, den gleichen 
Zustand cinzunchmcn. Zwei Fcrmionen mit der gleichen Energie miissen deshalb entgegengesetzten 
Spin haben. Auffallend ist, dass der Grundzustand"^ der Fermionen und Bosonen nicht bei E = 
liegt. 
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Temperatur besetzt werden. Diese Besetzungszahl n kann iiber das Grofikanonische Potential 
berechnet werden (Abb. 

Uber den Hamilton-Operator H fiir das aus iV-Teilchen bestehende Gas, kann die Dichtematrix 



B 




Abbildung 1.2: Das ideale Gas (A) befindet sich in einem viel groficren Warmebad (B), 
wodurch cine konstantc Temperatur erreicht wird. Ersteres ist damit im thermischen Gleichgewicht, 
wozu Energieaustausch vom Gas zum Warmebad notig ist. Die Pfeile zeigen an, dass zusatzlich 
Teilchenaustausch stattfindet. 



p des grofikanonischen Ensembles bestimmt werden 

1 



(3(H-fiN) , (1.1) 



p = —exp 

wobei Zq der Normierungsfaktor ist, der als Grofikanonische Zustandssumme bezeichnet 
wird. Im Exponenten treten noch der Teilchenzahloperator^ N, P = l/ksT mit der Temperatur 
T und das chemische Potential /.t in Erscheinung. 
Mittels Gl. ()1.2() ist es moglich, die mittlere Besetzungszahl n{e^ 

n{er) = r , / . (1.2) 

fiir einen Zustand mit der Energie ef bei einer bestimmten Temperatur T zu bestimmen. Oft 
wird in der Literatur eine GroBe z = exp{p/kBT) eingefiihrt, wodurch sich fiir die mittlere 
Besetzungszahl 

n{e^) = ie(^r) _ i (1.3) 

ergibt. 

Als Vorbereitung fiir den folgenden Abschnitt sei noch darauf hingewiesen, dass die Summation 



^Bosonen haben ganzzahligen Spin, hier werden jedoch nur Bosonen mit s = verwendet. 
^Fermionen haben halbzahligen Spin. 

^Der Grundzustand bezeichnet die niedrigst mogliche Energiekonfiguration des Systems. 
*Der Teilchenzahloperator angewendet auf einen Zustand liefert die Anzahl aller Teilchen. 
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der mittleren Besetzungszahl iiber alle moglichen Zustande die Teilchenzahl A'^ ergibt: 



N = ynr. (1.4) 



Es ist G1. H1.2|) zu entnehmen, dass mit abfallender Temperatur energetisch tiefer gelegene 
Niveaus bevorzugt werden. Um zeigen zu konnen, was passiert, wenn die Temperatur gegen 
Null geht, wird als Beispiel ein einfaches Potential verwendet. 

1.1 Das ideale Gas im harmonischen Potential 

Die Eigenwerte e^-* des Gases im dreidimensionalen harmonischen Potential ^muj'^x'^ sind 

ej-= IxJtUJx + ly^y + IzJ^i^z lx,y,z = 0, 1 ,2, 3. . . (1.5) 

wobei die Grundzustandsenergie eg = h/2 = {lOx + cOy + ujz) ist. 

Das weitere Vorgehen ist, e^^in Gl. (|1.4|) einzusetzen, mit dem Ziel, zu berechnen, wie sich 
die Teilchen bei der Annaherung an den absoluten Nullpunkt der Temperatur verhalten. 
Ein ausdruckstarkes Mafi dafiir ist das Verhaltnis von Teilchen im Grundzustand A'^o zu 
den im Gas enthaltenen Teilchen A^. Jene Teilchen, die nicht den Grundzustand bevolkern, 
bezeichnet man als anger egte Teilchen N' = N — Nq. 

Lasst man bei der Summation in der Formel fiir die Gesamtteilchenzahl (|1.4j) / = weg, 
ergibt sich fiir 

^' = E™r=E -exp{Plfuv)-l <Y,[expiPinw)-l]-'^N:^^, (1.6) 



1^0 1^0 



eine nach oben abgegrenzte Anzahl an angeregten Atomen N'^^^. Diese Naherung gilt, wenn 
1/ z gegen eins geht, also fiir T gegen Null. Zufiihrung von weiteren Teilchen ins Kondensat 
bewirkt eine direkte Anreicherung des Grundzustands. Die Grundzustandsenergie steckt im 
chemischen Potential ^. 

Nach einigen Schritten^ ergibt sich 

^0 - - " 



^Eine ausfiihrliche Herleitung findet sich in 
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mit dem als kritische Temperatur eingefiihrten Wert^ 

Unterhalb der kritischen Temperatur findet damit eine Teilchenanreicherung im Grundzustand 
statt. Im nachsten Abschnitt soli erklart werden, warum es sich dabei um einen Phaseniibergang 
handelt. 



1.2 Der Phaseniibergang und sein Ordnungsparameter 




A B C D 

Abbildung 1.3: Ob Teilchcn - odcr WcUcncharaktcr hangt von der Temperatur ab: Wahrend fiir 
hohe Temperaturen (A) die Atome als kleine Kugeln betrachtet werden konnen, nimmt das Atom fiir 
niedrige Temperaturen Wellencharakter an. Den Atomen kann eine thermische de Broglie Welle (B) 
zugeordnct werden. Dicsc Wellenlange nimmt mit abnchmcnder Temperatur zu. Bei Tc bcginncn die 
Wellenfunktionen sich zu iiberlappen (C). Angekommen beim absolutcn Nullpunkt der Temperatur 
bilden alle Teilchcn zusammen eine makroskopische Welle (D). 



Jedem Atom kann eine thermische de Broglie- Wellenlange At = \/ I-kT? j mkBT zugewiesen 
werden. Der Vergleich dieser Grofie mit den mittleren Abstanden / = ^^^^J^^ — im Gas zeigt, 
dass sich die de Broglie- Wellen oberhalb von nicht iiberschneiden, unterhalb schon. Bei der 
kritischen Temperatur selbst stimmen die thermische Wellenlange und der mittlere Abstand 
ungefahr iiberein. Dieses Anwachsen von Ay unterhalb bewirkt ein kollektives Verhalten 
der Atome im Gas. 

Ein Phaseniibergang ist dadurch charakterisiert, dass eine Grofie existiert, die unterhalb 
der kritischen Temperatur endlich ist und oberhalb der kritischen Temperatur verschwindet 

'^C(3) = 1.202... 
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|19j . Diese Grofie wird als Ordnungsparameter bezeichnet. Hier ist der Ordnungsparameter 
die makroskopische WellenfunktionJ 

Generell ist aber auch die Prage zu beantworten, ob die Annahme eines idealen Gases 
hinreichend genau ist. Dazu verwenden wir experimentelle Daten und vergleichen sie mit den 
mathematischen Vorhersagen dieses Abschnitts. 

1.3 Der Vergleich mit dem Experiment 

Wird der Bruchteil der Atome im Kondensat iiber die Temperatur aufgetragen und mit 
Gl. (|1.7j) verglichen, stimmt der theoretische mit dem gemessenen Wert hinreichend gut 
iiberein. Auch wenn bei der kritischen Temperatur aufFallt, dass T^^^ unterhalb liegt. 




Abbildung 1.4: Die Punkte stellen die Versuchsergebnisse von Ensher at al. |H] dar, wahrend die 
gcstrichclte Linic den matlicniatischen Vorhersagen fiir ein ideales Gas G1. H1.7|I entspricht. 



Dies aheine reicht jedoch nicht aus, um zu beantworten, ob die Wechselwirkung im Gas 
vernachlassigt werden kann. Ein probates Mittel zum Vergleich ist dagegen die Gesamt- 
energie E, da die Wechselwirkung Ejnt nur beriicksichtigt werden muss, wenn sie einen nicht 
verschwindenden Beitrag zur Energie E = Ekin + E^xt + Eint leistet. 

'^Der Exponent, welcher in Gl. 11.71 auftritt und als kntischer Exponent bezeichnet wird, weicht vom 
allgemeinen Fall No/N^tl- [T/T'J)'^^^ -unendlich grofier Potentialtopf - ab. 



7 



Im Versuch von Enscher et al. [S] wird ein Kondensat gebildet und das auBere Potential E^xt 
nach einiger Zeit abgeschaltet. Die reduzierte Energie E = Ekm + -E'mi wird in Abhangigkeit 
von der Temperatur gemessen. Die Durchfiihrung des Experiments ist jedoch mit der Schwierigkeit 
behaftet, dass sich ab dem Zeitpunkt der freien Expansion das Kondensat ausbreitet. Damit 
wachsen die Abstande zwischen den Atomen an, wodurch die Wechselwirkung abnimmt, 
deren Einfluss man messen mochte. 



CQ 



LU 




0.4 0.8 1.2 1.( 



1.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 
T / T^(N) 

Abbildung 1.5: Die gestriclielte Linie zeigt den theoretisch erwarteten Verlauf der Energieabnahme 
fiir ein ideales Gas. Dieser weicht von dem gemessenen Werten deutlich ab [H|. 



Oberhalb von ist eine gute Ubereinstimmung zu erkennen, das Gas kann als hinreichend 
ideal betrachtet werden. Die Situation andert sich aber fiir Temperaturen unterhalb von T^. 
Die Atome beginnen sich am untersten Energieniveau anzusammeln. Es wird also erwartet, 
dass sich die relative Energieabnahme 

NkgT NkgT \ N I keT \tJ kgT ^ 

verhalt. Betrachtet man jedoch Abb. H1.5() . ist eine lineare Abnahme zu erkennen. 
Damit muss die Annahme eines idealen Gases aufgegeben werden. Im nachsten Kapitel 
werden wir uns mit den nicht zu vernachlassigenden Wechselwirkungen beschaftigen. 
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Kapitel 2 

Das annahernd ideale Gas 



Die Widerspriiche zwischen Experiment und Theorie erfordern es, die Wechselwirkung zwischen 
den Teilchen des Gases zu beriicksichtigen. Es handelt sich dann um ein reales Gas. Im 
allgemeinen sind die Wechselwirkungen in einem N-Teilchensystem sehr kompliziert, da jedes 
Teilchen auf jedes andere Einfluss nimmt. 

Wir bewegen uns jedoch in Temperaturbereichen von T < l^K, in welchen fast alle sich im 
thermischen Gleichgewicht befindlichen Stoffe in der fliissigen Phase vorhegen. 

Dies kann nur umgangen werden, indem man zur Erzeugung von Bose-Einstein-Kondensaten 
extrem verdiinnte Gase {p < W^'^/cm'^) verwendet^. 

Die Wechselwirkung findet im verdiinnten Gas wesentlich reduziert statt. Es stofien immer 
nur zwei Teilchen gleichzeitig aneinander. Damit kann die Wechselwirkung zwischen den 
Atomen als kleine Storung zum idealen Gas betrachtet werden. Die Spezifizierung des Begriffs 
„verdiinntes Gas" und dessen mathematische Beschreibung ist Inhalt dieses Abschnitts. 



2.1 Ultrakalte verdiinnte Gase 

Es wird ein Gas betrachtet, dessen mittlerer Teilchenabstand groB gegeniiber der Wechselwirkungsreichweite^ 
ist. Damit konnen die Teilchen nur dann wechselwirken, wenn sie sich der Reichweite entsprechend 
annahern. 

^Zusatzlich muss die Kiihlung so rasch erfolgen, dass ein Erstarren verhindert wird. 

^Die Wechselwirkung zwischen Atomen kann proportional zu beschrieben werden, wobei r der 

Abstand sind. Damit ist die Reichweite unendlich. Aber es ist moglich zu zeigen, dass fiir Potentiale die 
starker als mit dem Abstand abfallen, die Annahme einer endlichen Wechselwirkung zulassig ist. 
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Abbildung 2.1: Das Bild zcigt cine schematischc Darstellung eines verdiinnten Gases: Die 
Atomabstande liegcn in einer GroBcnordnung von d ~ 1000 A^, die Wechselwirkungsrcichweite der 
Atome kann als endlich angcnommen werden und cs finden hauptsachlich Zwci-Tcilchcn-St5i3c statt. 



In die Rechnungen werden also nur Zwei-Teilchen-Wechselwirkungen einbezogen. 
2.1.1 Z wei-Teilchen- Wechselwirkung 

Wir betrachten zwei Atome im Schwerpunktsystem ()A.1|) . die im Abstand ri2 voneinander 
entfernt sind. Das Wechselwirkungspotential wird durch das Born-Oppenheimer- Potential 
fAbb. ()2.2|) ') beschrieben. In Abb. 1)2. 2(1 ist dargestellt, dass fiir niedrige Temperaturen gebundene 
Zustande moglich sind. Wiirde man das exakte Potential zur weiteren Berechnung verwenden, 
ware der Ubergang in die fliissige Phase unausweichlich. Es muss daher ein anderes Modell 
gefunden werden, das keine gebundenen Zustande erlaubt und eine weitere mathematische 
ErschlieBung ermdglicht^ . 

''W. Ketterle et aZ.|S| verwendeten zur Herstellung eines Bose- Einstein- Kondensats (1995) ein verdiinntes 
Sodium-Gas, bestehend aus bis zu 5- 10"" Atomen, mit einer Dichte von 10^'^ /cm^. Mit V ~ und p = No/N 
kann damit der mittlere Abstand berechnet werden. 

^Der exakte Potential- Verlauf ist fiir bestimmte Atome sehr schwer zu bestimmen und ein kleiner Fehler 
im Potential kann zu grofien Abweichungen im Resultat fiihren. 
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Abbildung 2.2: Dcr Potential- Vciiauf zwischcn zwci Atomcn wird vom rclativcn Abstand ri2 
bestimmt. Das Potential kann in drei Bcreiche eingetcilt werden. Fiir sehr kleine ri2 treten gebundene 
Zustande auf, ab rg streuen die Atome aneinander und im Abstand b voneinander doniiniert der Van- 
der-Waals-Anteil mit l/r^, endet die Reichweite der Wechselwirkung. 



2.1.2 Das exakte Potential 

Als Ausgangspunkt dienen die Ergebnisse der quantenmechanischen Streuprozesse, welche in 
Anhang HA.2|) kurz hergeleitet werden. 

9,, r pik\ri2-fi2'\ 

Airn'^ J \ri2 - ri2 '\ 

Diese Gleichung beinhaltet auch die unerwiinschten Bindungprozesse und Mehrfach- 
Streuprozesse. Eine Entwicklung von ^^(^2') = "00(^.2') + (^i^) bis zur Nullten Ordnung 
fiir schwache Potentiate schafft Abhilfe. Elastischen StoBe konnen liber 

Hri2) = e'^'^^ - / ^V{n2')e'^''' (2.1) 

47rn^ y \ri2-ri2'\ 

berechnet werden. Abb. (|2..S|) zeigt, dass die Streuung bei niedrigen Energien nur angewendet 
werden darf, wenn die Drehimpulsquantenzahl Null ist.^ 

^Die Drehimpulsquantenzahl zeigt an, wie die Elektronenschalen besetzt sind. Bei I — darf die letzte 
besetzte Schale kein p-Orbital sein. 
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Abbildung 2.3: Im allgemcincn Fall {I ^ 0) muss das V{r) (rot) auf cin cffcktivcs Potential Ve//(?') = 
—V{r) + Vi{r) erweitert werden. Der neue Term Vi = h^l{l + l)/2mr^ erlaubt es nicht mchr, dass fiir 
kleine positive Energien gebundene Zustande ausgeschlossen werden konnen. 



Beriicksichtigen wir, dass der mittlere Abstand ri2 im verdiinnten Gas grofi gegeniiber - 
der als endlich angenommen - Reichweite h ist, gilt ri2 <^ r^2- Folgende Naherung ermoglicht 
dies 

\n2-f'i2\=ri2-nf'i2 + o(J-y (2.2) 
wo n = die Streurichtung ist. 

ri2 " 

Fiir sehr niedrige Temperaturen liegen im thermischen Gas die auftretenden Energien in der 
Grofienordnung von ksT. Damit muss k in E = h^k'^/2ii klein sein, womit das Produkt 
/c6 <C 1 ebenfalls vernachlassigt werden kann. Unter diesen Voraussetzungen hangt der 
gestreute Zustand nicht mehr von der Streurichtung ab, well die Wellenlange groB gegeniiber 
den Details des Atompotentials ist. 

Ein Einsetzen der Annahmen - kleine Energien, endliche Reichweite im stark verdiinnten 
System - in Gl. (fTT|) . liefert 

V;(ri2) = e^'=^i2 (2.3) 
ri2 

Die neu eingefiihrte Variable a bezeichnet man als Streulange, welche mit der Streuamplitude 
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2l2 



I 
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(2.4) 



wie folgt zusammenhangt 



47r/i2 



2/i 



/ 



(2.5) 



Anzumerken ist nochmals, dass a nicht von der Streurichtung n abhangt, damit nur s- 
Wellen-Streuung betrachtet werden kann. 

Die hier beriicksichtigten Energien sind als sehr klein betrachtet worden. Lasst man 
nun k ^ gehen, nimmt die Energie des zu streuenden Zustands ab. Der Grenzfall eines 
Zustandes mit = 



kann zur numerischen Berechnung der Streuamplitude verwendet werden. 

Das exakte Potential kann mit Gl. (|2.6jl berechnet werden. Stellvertretend dafiir wird ein 
Modell- Potential eingefiihrt. Damit kann das Problem der gebundenen Zustande umgangen 
werden. 

Hinzu kommt, dass die Born- Approximation nur fiir schwache Potentiale giiltig ist, was fiir 
Abstande groBer als die Reichweite b des Potentials erfiillt ist, aber fiir kleine Abstande - 
dort wo gebundene Zustande auftreten konnen - nicht. 

Der Versuch, das starke Potential fiir sehr kleine Abstande durch hohere Ordnungen in der 
Born-Serie zu beseitigen, wiirde gebundene Zustande implizieren. 

Es wird ein Modell- Potential gesucht, das mit dem Verhalten - abgesehen von den Bindungszustanden 
- iibereinstimmt. 

2.1.3 Das Modell-Potential 

Die Idee fiir das Modell-Potential ist in Abb. ()2.4|) dargestellt. Es handelt sich dabei um ein 
Pseudo-Potential, 



1pE=Q = 1 

ri2 



(2.6) 



{rif2\Vint\'>Pi,2) = gS{ri - r2) 



d 



(n2V'i2(n,r2) 



(2.7) 



dri2 



J ri2=0 
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VAtompotentia 




Abbildung 2.4: Dargestellt ist ein beliebiges Potential, dessen Reichweite als endlich betrachtet 
werden kann. Dieses fiir kleine ri2 komplizierte Potential wird durch ein Kontaktpotential ersetzt, 
dessen Radius mit der Lage des letzten gebundenen Zustandes zusammenhangt. 



das sich aus einem Kontaktpotential und einem Regularisierungs-Operator zusammensetzt 
j^. Die Kopplungskonstante 

g = a (2.8) 

m 

beinhaltet die einzige Bestimmungsgroi3e - die Streulange - welche experimentell bestimmt 
wird. 

Die Streulange hangt nicht von der Reichweite des Potentials ab. Es besteht aber eine 
direkte Verbindung mit der Position des letzten angeregten gebundenen Zustands. Dessen 
Energieniveaus zweiatomiger Molekiile hangen von den Ubergangen zwischen den Schwingungszustanden 
des Molekiils ab. Dabei schwingen die beiden Atome um den metastabilen Abstand ro 
fAbb. ()2.2|) gegeneinander, es bildet sich ein sog. Vibrationspektrum aus. In der Vibarationsenergie 



E, = {{iyD-i^)K{ti,Ce)f/^ 



(2.9) 



kommen die Vibrations quantenzahl v und die beiden Konstanten K und vd vor. Die Vibrationsquantenzahl 
kann ganze positive Zahlen grofier Null = 0, 1, 2, .. annehmen und muss kleiner gleich V£, 
sein. Ein geradzahliges vy) entspricht damit dem letzten gebundenen Zustand, bei welchem 
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Abbildung 2.5: Dargestellt ist die Strculangc fiir Rb2 (im Triplctt-Zustand) Potential in 
Abhangigkcit von vo- Die Divergenzen stellen sich fiir geradzahligc h' = cin. Grofic negative Werte 
von a treten dann auf, wenn das Potential den nachsten gebundenen Zustand gerade vermeiden kann. 
Fiir grofie positive Werte ist es genau umgekehrt. Sie liegen vor, wenn das Potential stark genug ist, 
den letzten gebundenen Zustand zu unterstiitzen |15j . 



gleichzeitig die Trennung des Molekiils eintritt |15j . 
Eine exakte Berechnung der Streulange a ergibt 



mit der Konstanten 



und 



:a,,(l-tan(<^)), (2.10) 
r(3/4) /2/xC6 



2V2r(5/4) V 



1/4 

(2.11) 



(/. = 7r(i/D + 1/2). (2.12) 

Mit und der Atommasse ist die Streulange fixiert. Durch Manipulationen an der Tiefe 
des Potentials und der Masse ist es moglich, die Streulange zu variieren (Abb. H2.5() ). 

In alien weiteren Berechnungen wird als Wechselwirkungspotential dieses Modell-Potential, 
ein Kontakt-Potenial mit der Streulange a als Radius, verwendet. Damit ist es moglich, die 
lokal abhangigen Wechselwirkungen durch ein mittleres Feld zu ersetzten. 
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Kapitel 3 

Die Gross-Pitaevskii-Gleichung 



3.1 Der Viel-Teilchen-Hamilton-Operator 

Die quantenmechanische Beschreibung eines Viel-Teilchen-System - in dem nur paarweise 
Wechselwirkung stattfindet - erfolgt im Formalismus der zweiten Quantisierung^ . In diesem 
ist es moglich, Teilchen zu erzeugen bzw. zu vernichten. Einer Anregung aus dem Grundzustand 
entspricht zum Beispiel eine Vernichtung eines Teilchens im Grundzustand, kombiniert mit 
einer Erzeugung eines Teilchens im angeregten Zustand. Diese Beschreibung eignet sich zur 
Analyse eines Anregungsspektrums. 

Das Verhalten eines Teilchens an einem bestimmten Ort liefert der sog. Feldoperator ^^(x), 
der ein Teilchen bei x erzeugen bzw. vernichten kann. 

Als Ausgangspunkt verwenden wir den groiBkanonischen Hamilton-Operator K = H — 
im Feldformalismus 



Aui3er der Zwei- Teilchen- Wechselwirkung wurden bislang keine weiteren Annahmen gemacht. 
In den nachsten Schritten werden zusatzliche Naherungen fiir das ultrakalte verdiinnte Gas 
durchgefiihrt. Das chemische Potential ist in Vext enthalten. 

^Zweite Quantisierung, weil analog zur ersten Quantisierung, bei welcher die Observablen durch Operatoren 
ersetzt wurden, nun auch noch die Zustande durch Operatoren beschrieben werden. 




X - x')^{t,x')^{t,x). 



(3.1) 
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3.2 Der Bogoliubov-Ansatz 



Hierbei handelt es sich um ein Verfahren zur Berechnung des Anregungsspektrums fiir 
Temperaturen T <^ Tc- Es wird davon ausgegangen, dass sich fast alle Atome im Grundzustand 
befinden, also N — Nq <^ Nq ~ N. Der Grundzustand ist die makroskopische Wellenfunktion, 
die sich unterhalb der kritischen Temperatur ausbildet und durch die komplexwertige Funktion 

in der Naherung fGl. (|3.2() ) zur Geltung kommt. Anregungen konnen dann als Fluktuationen 
um den Grundzustand aufgefasst werden und sind durch den Feldoperator 6^{x) beschrieben. 



Das e hebt hervor, dass die Anregungen klein sind. 

Bevor nun mit diesem Ansatz weiter gerechnet wird, bedarf es einer Klarung beziigUch 
der komplexwertigen Funktion 

3.2.1 Spontane Symmetriebrechung 

Wie schon erwahnt, beschreibt ^(x) den Kondensatanteil. Es wird angenommen, dass der 
Feldoperator einen nichtverschwindenden endlichen Erwartungswert (^) = ^(x) hat. Der 
Erwartungswert der Fluktuationen sollte damit verschwinden {6^) = 0. Die Annahme eines 
nichtverschwindenden Erwartungswertes bringt jedoch eine Symmetriebrechung mit sich. 
Es handelt sich um die Eichsymmetrie, welche verletzt wird. Diese fordert, dass der Hamilton- 
Operator invariant gegeniiber derartigen 



Transformationen mit der Phase Oq ist. 

Die Forderung, dass der Erwartungswert (0|^'(x)|0) des Grundzustand |0) nicht verschwinden 
darf verletzt die Eichinvarianz und damit die Symmetric. Konvention ist, a = zu wahlen. 
Bei der Behandlung der idealen Gase wurde gezeigt, dass es sich um einen Phaseniibergang 
handelt. Da jeder Phaseniibergang zweiter Ordnung mit gebrochener Symmetric verbunden 
wird, muss es sich spontan um einen solchen handeln. 



^'(t, x) = ^{t, x) + eil){t, x) + ... . 



(3.2) 



^{x) ^{x)e' 



(3.3) 
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3.3 Die Mean-field-Naherung 

Nachdem die Bogoliubov-Naherung diskutiert wurde, kann diese in den Hamilton-Operator 
fGl. (|3.1|l ) eingesetzt werden, wobei zunachst fiir das Wechselwirkungspotential das Kontaktpotential 
eingesetzt wurde V{x — x') = g6{x — x'). Dabei ergeben sich folgende Terme 

k = J <fx{<^*{x)kQ^{x) + f 

+e j d^x (^*{x)KQ5^{t,x) + g^*{x)^*{x)^{x)5'^{t,x) + konjungiert komplex^ 

+£2 / d^x x)koS^{t, x) + f ^>*(x)2<5^'(t, xf + ^<^{xf6i'^{t, xf + 2g^xf6i/Ht, xf^ 

+e^g f d^x (^^*{x)6'if^t,x)S'^{t,x)6'if{t,x) + $(f)J^'t(t, f)5^t(t^ f)j 

+e^f / d^x (^6if^{t,x)64!^{t,x)64/{t,x)64!{t,x)^ 

(3.4) 

wobei ko = (-^V2 + Ve^tix) - /u) ist. 

Die Terme dritter und vierter Ordnung in e konnen durch eine Mean-Feld-Naherung 
ersetzt werden. Dazu wird das lokale Feld durch seinen Mittelwert, also durch ein mittleres 
Feld ersetzt. Zum Beispiel wird aus 6^^t,x)6^{t,x)6^{t,x) 2(6^'^ {t, x)5^ {t, x))6^ {t, x) + 
(5^(t, x)(5^(i, x))5^'^(t, x). Diese Naherung wird als Hartree-Fock-Bogoliubov-Naherung (HFB) 
bezeichnet [T3] . 

Vernachlassigt man weiterhin anomale thermische Dichten (6^6^), bekommt man die Popvo- 
Naherung, die fiir Temperaturen bis zu 0.7 Tc verwendet werden kann. Vernachlassigt man 
weiterhin alle thermischen Dichten ((5^5^^) erhalt man die Gross-Pitaevskii-Gleichung, die 
als nachstes untersucht wird. 



3.3.1 Die zeitunabhangige Gross-Pitaevskii-Gleichung 

In diesem Abschnitt wird der Grundzustand des fast idealen Gases mit dem Ritzschen 
Variationsprinzip berechnet. Die Forderung, dass alle Atome das niedrigste Energieniveau 
besetzen, wird als Nebenbedingung 



d^x\^{x)\^ = 1 (3.5) 
mit dem Lagrange Parameter A im Energiefunktional 

(if^fe) = {H) - XNo (^j d3x|^>(x)|2 - 1^ (3.6) 
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beriicksichtigt. 

Als Hamilton-Operator kann Gl. ()3.4|) verwendet werden, wenn man bedenkt, dass sich der 
groBkanonische Hamilton-Operator zu diesem reduziert, sobald alle Glieder mit fi weggelassen 
werden. Mit (5^) = fallen bei der Berechnung des Erwartungswerts alle Terme mit e weg^. 
Das sich ergebende Funktional 

^[$,$*]= I d^x (^<^*{x) (^-^V^ + Ve.t{x)^ c|,(f) + ||$|4^ -AiVo^y d3x|$(x)|2- 1^(3.7) 

wird um den Grundzustand <I>(x) ^{x)+6^{x) variiert. Fordert man, dass die Funktionalableitung 
6/6^ der Grofie 5E = E[^ + 6^, <I>* + 5$*] fiir alle 6^ verschwindet, muss folgende Gleichung 
erfiillt sein: 

r ?,2 ^ .1 

<^{x) = 0. (3.8) 



-7^V2 + Vecctix) + iVo5|^(x)|2 - A 

Dieser Ausdruck ist - mit der physikalischen Interpretation des Lagrange-Multiplikators als 
chemisches Potential A = /U - als zeitunabhangige Gross-Pitaevskii-Gleichung bekannt. Als 
Losung ergibt sich der Grundzustand ^>(x) = ^>o(^^)- 

Die Frage, ob diese Naherung die Realitat besser beschreibt als die theoretischen Vorhersagen 
fiir das ideale Gas, klart Abb. 1)3. 1(1 . Sie zeigt den Vergleich zwischen idealem Gas und dem 
ultrakalten verdiinnten Gas (Beschreibung nacli Gross-Pitaevskii) mit dem Experiment. 

3.3.2 Die zeitabhangige Gross-Pitaevskii-Gleichung 

Bei der Herleitung fiir den zeitabhangigen Fall wird nicht auf die Variationsrechnung zuriickgegriffen 

3 

Der Hamilton-Operator inklusive Zwei-Teilchen-Wechselwirkung in Feldquantisierung ist 

H= dx^\t,x)[-—V^ + Ve^t(.x)]4f{t,x) + 

dxdx' i/{t,x)^{t,x')V{x - x')^\t,x')¥{t,x). (3.9) 

Die Zeitentwicklung eines Feldoperators, auf den H wirkt, kann direkt durch die Heisenberg- 
Relation 

ih^^^{t,x) = [4>{t,x),H] (3.10) 



^Die Terme dritter und vierter Ordnung fallen erst nach der Mean-Feld-Naherung weg, weshalb Gl. 13.41 
nach der Einfiihrung des Mittleren-Feldes zu verwenden ist. 

^Die Herleitung der zeitabhangigen Gross-Pitaevskii-Gleichung mittels des Variationsprinzips findet sich 
inm. 
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-60 -40 -20 20 40 60 

z(/im) 

Abbildung 3.1: Die expcrimcntcllcn Datcn (Punktc) stimmtcn nicht mit den thcorctischen 
Berechnungen iiberein, falls das Gas als ideal angcnommcn wird (gcstrichcltc Linie). Die 
durchgezogene Linie zeigt die thcoretische Vorhersagc dcr Gross-Pitacvskii-Glcichung )14| . 



angegeben werden. Setzt man im Anschluss daran fiir den Feldoperator ein klassisches Feld 

^ ^ \/iVo^ (3.11) 
fiir den Grundzustand ein, ergibt sich die zeitabhangige Gross-Pitaevskii-Gleichung: 



2m 



(3.12) 
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Kapitel 4 

Anregungen im Kondensat 



In diesem Abschnitt werden kleine Anregugen im Bose-Einstein-Kondensat untersucht. 
Diese konnen als Schallwellen interpretiert werden, deren Ausbreitungsgeschwindigkeit als 
Schallgeschwindigkeit bezeichnet wird. 

Der erste experimentielle Nachweis von Schallwellen im Bose-Einstein-Kondensat konnte 
erstmals in einer frei expandierenden Atomwolke erbracht werden. Dieser Versuch wird 
am Ende des Abschnitt kurz vorgestellt. 

4.1 Das Anregungsspektrum 

Fiir T = konnen alle auftretenden Dichten Gl. (|3.4|) vernachlassigt werden.^ Nimmt man 
zusatzlich an, dass sich die gesamte Dichtemodulation mit einem festen Impuls p = hVO im 
Kondensat ausbreitet (siehe Abb. (|4.1|) ). ergibt eine langere Rechnung das folgende Anregungsspektrum: 




(4.1) 



Eine Herleitung fiir Gl. ()4.1|) ist in zu finden. Ftir dieses Energiespektrum muss die 
Dichteschwankung lokal auf einen sehr kleinen Raum begrenzt sein, was fiir Fluktuationen 
kleiner als die charakteristische Lange des Oszillators a_f/o = (h/mio)^^'^ erfiillt ist. Zusammen 



^Unter dem BegrifT einer kleinen Storung versteht man eine lokale Dichtemodulation, welche den BEC- 
Zustand nicht vernichtet und als koUektive Anregungen aufgefasst werden. 
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A) quantenmechanisches Bild 




T 



B) klassisches Bild 




Abbildung 4.1: In der Quant enmechanik (A) kann sich keine Dichteverteilung fortbewegen, ohne 
deformiert zu werden. In der klassischen Mechanik (B) ist es dagegen moglich, die Storung als ganzes 
mit einem Impuls propagieren zu lassen. 



mit der Unscharferelation ergibt ApAx > ^ fiir den Impuls der Welle p ^ h/auo- Der 
Radius Rq ~ 2kBT /muP' einer thermischen Atomwolke hangt von der Temperatur ab. 
Zusatzlich muss noch fiir den Impuls Rq » ano gelten, und damit T S> hio/kB- Zu 
Anfang des Abschnitts wurde eine Temperatur nahezu T <^Tq vorausgesetzt. Durch hohere 
Teilchenzahlen wachst Tq an, wodurch beide Bedingungen erfiillt werden konnen. Im Thomas- 
Fermi-Limes Na/ttHO ist die kinetische Energie im Vergleich zur Wechselwirkungsenergie 
innerhalb des Kondensats sehr klein. Ohne den kinetischen Term T4xt— ^ = gpo{^) vereinfacht 
sich das Energiespektrum Gl. (|4.1|) zu 

E'^{p,x) = ^{p' + 4mgpo{x)f^ . (4.2) 
2m 

In Abb. (|4.2|) ist das Energiespektrum in Abhangigkeit vom Impuls dargestellt. 



4.2 Kleine Anregungen im Kondensat 

Fiir kleine p-Werte lasst sich die Dispersionsrelation einfach angeben E = pc wobei 

V m 

die Schallgeschwindigkeit ist. 
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Impuls p 



Abbildung 4.2: Die blauc Kurve stellt den Energicverlauf in Abhangigkeit vom Impuls dar. Die gelbe 
Gerade ist zur Orientierung eingezeichnet und zeigt, dass in zwei Bereiche einzuteilen ist. Fiir kleine 
Impulse ergibt sich ein lineare. damit ein phononenartiges Spcktrum mit E = pc. Dcr quadratische 
Ansticg fiir grofiere Impulswerte lasst auf Teilchenanregung E = /2m + Eq schlieficn. 

4.3 Experimentelle Beschreibung 

Ketterle et al. [2] haben die Ausbreitung von Schall im Bose-Einsteinkondensat gemessen. 
Im Versuch wurden 5 x 10^ Natrium- Atome abgekiihlt und durch ein geeignetes Magnetfeld 
wurde die Atomwolke zigarrenformig angeordnet. Die Dichtestorung wurde im Zentrum 
des Kondensats durch einen Laser erzeugt. Die Anregung propagierte dann spharisch vom 
Zentrum nach aufien, was direkt gemessen werden konnte. 
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Abbildung 4.3: Ausbrcitung einer kleinen St5rung im Bose-Einstcin-Kondcnsat. Zunachst wird das 
Kondcnsat in cincr magnctischen Falle erzeugt. Zur Zeit t ~ wird cin Laser kurzzeitig an (a)- 
bzw. ausgcschaltct (b). Damit wird cine positive/ negative Dichtemodulation erzeugt, die mit der 
Schallgeschwindigkeit durch das Kondensat propagiert 
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Abbildung 4.4: Bcobachtung dcr Schallwellen in cincm Kondcnsat. Zur Mcssung wurdc cine das 
Kondcnsat nicht zcrstorcndc Phascn-Kontrast- Abbildung verwendet. Die cinzelncn Plots zeigcn das 
Kondensat im Abstand von 1.3 fis. In bcidcn Vcrsuchcn wurdc das Kondcnsat mit cincm Lascrpuls 
- nach Anregung der Storung - in zwei Telle getrennt. Die Abbildungen zeigen zwei verschieden 
praparierte Kondensate: Die Lange der oberen Kondensatwolke ist 450 fj,m. In der unteren Sequenz 
ist der radikale Einschluss geringer als im Kondensat zuvor, was zu cincm radialen Ausweiten der 
Dichte fiihrt. Es ist dcutlich zu schen, dass die Schallgeschwindigkeit im untcrcn Bild langsamer ist 
als im oberen. Dort ist die Dichte, wclche proportional zur Schallgeschwindigkeit ist, geringer [2]. 
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Teil II 



Einfiihrung in die Allgemeine 
Relativitatstheorie 
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Kapitel 5 

Eigenschaften schwarzer Locher 



5.1 Allgemeine Relativitatstheorie 

Die Allgemeine Relativitatstheorie beschreibt gravitative Felder. Die Dynamik ergibt sich 
aus der Einsteinschen Feldgleichung 

= -kT-^P. (5.1) 
Rechts ist der Energie-Impuls-Tensor der Materie 

mit der Determinante der Metrik g = det{gap) und der Materie- Lagrangefunktion Lm = 
LM{9ai3\d'ygaf3\^^\^'te) abgebildet. fasst alle Materievariablen zusammen. Auf der linken 
Seite von Gl- HS-lf) steht der Einstein- Tens or 

Gap = Rafi — -^Rgaf^, (5-3) 

mit dem Riemannschen Kriimmungstensor Ra/s = '^{gs-y, 

und dem Ricci-Skalar R = Rapg"^^ ■ Die Konstante k = beinhaltet die universielle 
Gravitationskonstante G ist [T^ . 

Die Einstein- Gleichung ist eine Bestimmungsgleichung fiir die Metrik. Die Vorgehensweise 
ist, zunachst die Materieverteilung zu bestimmen, in Gl. (|5.1() einzusetzen und mit der zweiten 
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Einstein- Gleichung^ 



Tal3;l3 = 



(5.4) 



die Metrik zu ermitteln [Jj. 

Einer bestimmten Massenverteilung wird eine Metrik 



9al3 



900 901 902 903 

901 9n 912 913 

902 912 922 923 
\ 903 913 923 933 ) 



(5.5) 



zugeordnet, die sich aus den Einsteinschen Feldgleichungen ergibt. Die Eintrage in der Metrik 
Qa^ji = ga/3{x) sind im Allgemeinen von der Zeit und dem Ort abhangig^. 

Die Metrik ist ein symmetrischer Tensor zweiter Stufe, der es ermoglicht, den infinitesimalen 
Abstand ds zweier benachbarter Punkte x und x + dx im gekriimmten Raum zu bestimmen: 



ds' 



9, 



^I3{x)dx°'dx^ . 



(5.6) 



Fiir den flachen Raum sind die Eintrage in Gl. (|5.5() unabhangig von der Raum-Zeit. In 
unserem Universum kann der flache Raum durch die Minkowski- Metrik 



/ 1 \ 

0-10 

0-10 

-1 



(5.7) 



beschrieben werden. Das Linienelement fiir den flachen Raum ist 



ds' = c'{dx'>f - {dx'f - {dx'f - {dx') 



1\2 



„2\2 



(5. 



In dieser Arbeit sind zwei bestimmte Metriken von Interesse, welche im Folgenden vorgestellt 
werden. 



'^Wobei das Semikolon die kovariante Ableitung nach (5 bezeichnet. 

■^Raum und Zeit sind in der Relativitatstheorie zur Einheit verschmolzen und werden zu einem Vierervektor 
zusammengefasst. Die nullte Komponente des Vierervektors ist die Zeit, die drei verbleibenden Eintrage 
entsprechen den dreidimensionalen Ortskoordinaten. 
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5.1.1 Die Metrik fiir ein schwarzes Loch 



Die exakte Losung eines kugelsymmetrischen, statischen Himmelskorpers wurde zuerst von 
Schwarzschild berechnet. Die Dynamik fiir Teilchen aui3erhalb einer solchen Materieverteilung 
ist in spharischen Koordinaten (t, r, 6, (p) gegeben durch 

„2 _ 2 A 2m \ 2 dr"^ „2 .o2 



ds' = c'[l-—Jdt'- - r'dn', (5.9) 

mit dem Raumwinkelelement dO, = sin[9)d9dip und m = GM/c^, wobei M die Gesamtmasse 
der Materieverteilung ist. 

In der Schwarzschild-Metrik treten zwei Besonderheiten auf. Fiir r = 2?7i = rg verschwindet 
gooii^En) = 0. rEH wird als Gravitationsradius bzw. Schwarzschildradius bezeichnet. Innerhalb 
r < rg^ ist §00 < und damit der Tangentenvektor - Xjcdt' an die Weltlinie ds^ = 
(l — ^) dt"^- raumartig, fiir goo{fg) = liclitartig und ausserhalb r > Vg ist 500 > 
zeitartig. 

Allgemein wird eine Stelle im Raum, an welcher der Tangentenvektor sein Vorzeichen andert 
als Ereignishorizont bezeichnet. Der Ereignishorizont ist eine zeitartige Hyperflache, d.h. der 
Tangentenvektor am Ereignishorizont ist parallel zur Zeit-Achse. In Schwarzschild-Koordinaten 
wird gii = , \,n\ am Ereignishorizont unendlich. Ein sich dem Ereignishorizont radial 

V r ) 

nahernder Korper mit der Geschwindigkeit 



hm — = — 

TEH 



wird auf f = abgebremst. Damit ist es unmoglich, in ein schwarzes Loch einzutreten. 
Eigentlich wird erwartet, dass durch die mit abnehmender Entfernung steigende gravitative 
Anziehungskraft die Geschwindigkeit zu- und nicht abnimmt. Bei der physikalischen Interpretation 
muss die Wahl des Koordinatensystems beriicksichtigt werden. Die Schwarzschild-Koordinaten 
entsprechen einem Bezugssystem aufierhalb des schwarzen Lochs. Von dort wird ein Korper - 
z.B. ein Raumschiff - beobachtet, wie dieses sich dem schwarzen Loch nahert. Ein Beobachter 
am Ereignishorizont - mit einem anderen Bezugssystem - sieht das Raumschiff in endlicher 
Zeit den Horizont passieren Bei r = rEH handelt es sich unter diesem Gesichtspunkt um 
keine Singularitat im eigentlichen Sinn, da sie mit einer geeigneten Koordinatentransformation 



^Das Linienelement Gl. 15.91 ist nur fiir den Bereich teh < r < 00 giiltig. 
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zu beseitigen ist. Am Punkt r = hingegen liegt eine echte (intrinsische) Singularitat vor, 
die in alien Bezugssystemen zu finden ist. 



5.1.2 Die Metrik fiir das de-Sitter-Universum 



Das de- Sitter- Modell beschreibt unser Universum approximativ, indem die Gesamtmasse des 
Universums gegenueber dessen Grosse verschwindend betrachtet wird. Das Universum kann 
daher als masselos und global flach betrachtet werde. Es kann berechnet werden, wie sich ein 
solches Universum mit der Zeit verandert. Eine Herleitung findet sich in [7]. Es ergibt sich 
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exp{2Ht) 



















(5.10) 



exp(2iJt) J 



die als de- Sitter- Metrik bezeichnet wird. Hier ist 



H= - 



1/2 



(5.11) 



mit der kosmologischen Konstanten A. Der allgemeine zeitabhangige Fall wird durch das 
Hubble- Gesetz 

m 



H{t) 



R{ty 



(5.12) 



beschrieben, wobei H{t) als Hubbel-Konstante bezeichntet wird. R{t) ist ein Skalenfaktor. 
Multipliziert mit einem, zum Zeitpunkt t = gemessenen Abstand, gibt er an, wie sich dieser 
mit der Zeit andert. 

Im de-Sitter-Modell ist A > 0. Es erfahrt eine exponentielle Expansion. ^ 



5.2 Feldquantisierung im Gravitationsfeld 

In diesem Abschnitt wird kurz erklart, was unter dem Begriff Hawking- Strahlung verstanden 

wird. Eine ausfiihrliche Diskussion findet sich in 

*Es wird noch darauf hingewiesen, dass es sich hier um das einfachste Modell fiir das Universum handelt. 
Auch wenn die Dichte der Materie im Universum zum jetzigen Zeitpunkt verschwindend ist, lasst sich diese 
Aussage nicht mehr aufrecht erhalten, wenn man lange genug in der Zeit riickwarts wandert. Auch die 
Hintergrundstrahlung, die im de-Stitter-Modell berechnet werden kann, entspricht nicht dem experimentell 
gemessenen Wert. 
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Abbildung 5.1: Die Entwicklung nach den Moden entspricht einem Wechsel vom Ortsraum (oberes 
Bild) in den Impulsraum (unteres Bild). Im Impulsraum sind die Moden mit k und —k miteinander 
gekoppelt, so dass nur paarweise Anregung von Moden mSglich ist. 



Bei der Quantisierung der Bewegungsgleichung eines massebehafteten Teilchens im gekriimmten 
Raum nimmt der Hamilton-Operator die Gestalt unendlich vieler entkoppelter Oszillatoren 
an. An jedem Punkt des Raumes kann man sich einen harmonischen Oszillator im Grundzustand 
vorstellen. Ein Teilchen am Ort x zu erzeugen, heifit den Grundzustand auf das nachstgelegene 
Energieniveau anzuregen. Ein Wechsel in den Impulsraum mittels Fourier- Transformation 
liefert einen Hamilton-Operator, welcher durch die sog. Bogoliubov- Transformation diagonalisierbar 
ist. Es ergibt sich, dass nie ein einzelnes Teilchen angeregt werden kann. Es tritt Paarteilchenerzeugung 
auf, wobei diese einen gleichgrofi entgegengesetzten Impuls haben. Mit anderen Worten, der 
Gesamtimpuls des Systems bleibt bei der Teilchenerzeugung erhalten (siehe Abb. 1)5. 1(1 ). Ob 
Teilchen erzeugt werden oder nicht hangt vom Gravitationsfeld ab. 

5.2.1 Hawking-Strahlung 

Unter Hawking-Strahlung versteht man Paarteilchenerzeugung am Ereignishorizont eines 
schwarzen Lochs. Es sind zwei physikalische Phanomene, die dafiir verantwortlich sind. 
Zum einen sagt die Unscharferelation AE-At > h voraus, dass iiberall im Raum Vakuumfluktuationen 
auftreten. Darunter versteht man, dass standig Teilchenpaare erzeugt werden, die sich nach 
einer gewissen Zeit At wieder vernichten. Es sei denn, es existiert eine aussere Kraft, die stark 
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genug ist, diese Teilchen voneinander zu trennen. Am Ereignishorizont ist diese Bedingung 
erfiillt. Ein Teilchen (mit Impuls k) wird als Hawking-Strahlung ins Weltall gesendet, das 
andere (mit Impuls —k) wird in das schwarze Loch gezogen. Dadurch verliert das schwarze 
Loch Energie, so dass die Teilchenerzeugung nicht ohne Aufwand betrieben werden kann. 

Eine wichtige Eigenschaft Schwarzer Locher ist, dass sie dynamisch stabil (siehe Abschnitt 
(|7.2)l ) sind. Das bedeutet, dass wenn man bei den Betrachtungen komplexe Frequenzen 
zulasst, die Zeitskala fiir das dynamische Anwachsen der dadurch produzierten Teilchen groi3 
ist. Schwarze Locher werden bisher als weitgehend stabil betracht. 
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Teil III 

Das Bose-Gas als 
Gravitationsmodell 
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Kapitel 6 



Die Wellengleichung fur die 
Storung 



Die Verbindung der Allgemeine Relativitatstheorie mit dem Bose-Einstein-Kondensat wird 
in diesem und im nachsten Kapitel erlautert. 

Im Folgenden wird fiir eine im Kondensat propagierende Storung die Wellengleichung aufgestellt 
(siehe Kapitel @) und gezeigt, dass es sich dabei um eine Wellengleichung in einem effektiv^ 
gekriimmten Raum handelt. 

Die Ergebnisse wurden bereits von Cirac et al. ^7] publiziert. Es wird aber darauf hingewiesen, 

dass hier der allgemeine Fall, dass die Kondensat grossen c = c{t, x) und v = v{t,x) vom Ort 

und der Zeit abhangen. In der zitierten Veroffentlichung waren die Schall - und Hintergrundgeschwindigkeit 

nur vom Ort abhangig. ^ 



^Der Begriff effektiv weist darauf hin, dass es sich nicht im eigentlichen Sinne um einen gekriimmten Raum 
handelt, weil dieser die Anwesenheit eines echten Gravitationsfeldes fordern wiirde. Trotzdem lasst sich eine 
Analogie finden, da in den Gleichungen ein symmetrischer Tensor zweiter Stufe auftaucht, der als effektive 
Metrik aufgefasst werden kann. Die ganze Information des gekriimmten Raumes ist in der Metrik enthalten 
und damit ist es zulassig, von einem effektiv gekriimmten Raum zu sprechen. 

^In Kapitel llUl wird ein zusatzlich von der Zeit abhangiges System verwendet. 



39 



6.1 Die hydrodynamischen Gleichungen 

Ausgehend von der zeitunabhangigen Gross-Pitaevskii-Gl. ()3 . 1 2 (1 wird eine kleine Variation 
in der Dichte p = po + £Pi und in der Phase 6 = 6q + eOi ^ in den Zustand <^ = y/pe^^ 
eingesetzt 

V 2m J (6.1) 

wobei pi und 9i die Storung beschreiben und das chemische Potential p in Vext enthalten 
ist. Im Anhang ist Gl. H6.ip ausgerechnet worden. Es ergaben sich fiir die Storung: 

(|n3|:0i=-A(v^oV^i)-^Pi 

m n 

^ : pi = - -V(poV^i + piV^o) 
m 

6.1.1 Die Wellengleichung 
G1. HB.3|) wird nach pi aufgelost 

Nogpi = -{vVei)-9i 
und anschliei3end nach der Zeit abgeleitet 

Nogpi = -dt{vV9i)-6i. 

Eingesetzt in G1. ()B.5|1 hefert 

-01 - dt{vV9i) - V{v9i) + V ((c^ - v'^)V9i) = 0, (6.2) 
eine Differentialgleichung fiir die Phase der Storung. 



Dass es sich dabei um eine Wellengleichung handelt, wird sichtbar, wenn 



\ 



-(? + vl 



Vz 



\ 



(6.3) 



Das e zeigt die Ordnung an, wir betrachten damit eine Entwicklung bis zu ersten Ordnung. 
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eingefiihrt wird. Jetzt lasst sich Gl. ()6.2|) schreiben als: 



^>.{^/^g^''^,e^) = (6.4) 

Formal handelt es sich um eine klassische Feldgleichung fiir ein masseloses Teilchen 9i in 
einem gekriimmten Raum. Die Rolle des gekriimmten Raumes nimmt dabei das Kondensat 
ein, weil samtliche Eintrage in der Metrik aus den GroBen des Kondensats bestehen. Es tritt 
die Schallgeschwindigkeit c fiir die Dichtefluktuationen aus dem vorhergehenden Abschnitt 
sowie die stationare Hintergrundgeschwindigkeit v auf. 



Fiir ein verschwindendes v, wird Gl. ()6.3() 





( 


1 











\ 


-1 

















^~ 






























-c2 


/ 



(6.5) 



zur Minkowski-Metrik, welche den flachen Raum beschreibt. Die Minkowski-Metrik fiir den 
flachen Raum, eingesetzt in Gl. (|6.4|) . liefert die gewohnliche Wellengleichung: 

{c'd^ - v2) 01 = 



6.2 Das efFektive Schwarze Loch 



Mit g^vQ^'^ = 1 kann die kovariante Metrik zu Gl. H6.3p berechnet werden. Es ergibt sich: 
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1 



(6.6) 



Es ist leicht zu liberpriifen, ob in einem System mit dieser effektiven Metrik ein schwarzes 
Loch auftritt, oder nicht. Aus Kapitel Q ist bekannt, dass ein Ereignishorizont vorliegt, wenn 
an einem Punkt goo = wird. Hier ist goo = —c^+v, damit ist es vom System abhangig, ob ein 
Ereignishorizont auftritt. Der Bereich, in dem die Hintergrundgeschwindkeit schneller fliefit 
als die Schallwellen sich ausbreiten konnen, kann mit dem Inneren eines schwarzen Lochs 
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in Verbindung gebracht werden. Dort konnen die Schallwellen nur mit der Richtung des 
Kondensats propagieren. Aui3erhalb, im Bereich wo die Hintergrundgeschwidikgeit kleiner 
als die Schallgeschwindigkeit ist, kann sich der Schall in beide Richtungen bewegen. 
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Kapitel 7 

Die Fluktuationen um den 
Grundzustand 



Die Untersuchung der Zeitentwicklung von Fluktuationen um den Grundzustand ist Bestandteil 
dieses Kapitels. 



7.1 Der Bogoliubov-Ansatz bis zur ersten Ordnung 

Der Ansatz von Bogoliubov ^{t,x) = ^{t,x) + e'4>{t,x) + ... beschreibt ein System, das 
sich bis auf kleine Fluktuationen im Grundzustand befindet. Bisher wurde die erste Ordnung 
nicht beriicksichtigt. Jetzt ist gerade diese von Interesse. Dazu greifen wir die Rechnungen aus 
(|3.3|1 auf. Anstatt der Mean-Field-Naherung wird als nachstes K in die Heisenberg-Relation 

ih^^if{t,x) = [i'{t,x),k], 

eingesetzt. Beriicksichtigt man zusatzlich die Vertauschungsrelationen der Feldoperatoren fiir 
Bosoen 



i^ix),ip{x') 



0, 



V^t(x),'0t(f' ) 



0, 



i^{x),ip\x')\ =6^^\x-x') (7.1) 



ergibt sich in erster Ordnung von e die zeitabhangige Gross-Pitaevskii-Gl. (|3. 12|) . Diese verschwindet, 
da Fluktuationen um den Grundzustand betrachtet werden und damit die klassische Wellenfunktion 
im Bogoliubov-Ansatz gerade dem Grundzustand der der makroskopischen Wellenfunktion 
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entspricht. Die Terme dritter und vierter Ordnung konnen vernachlassigt werden, so dass nur 
Ausdriicke mit zuriick bleiben. ^ Die zeitabhangige Bewegungsgleichung fiir die Fluktuationen 
ist 



•f,^ / 



2m 



+ Ve^t + 2Ngps 



■0 + Ng^'^ip^. 



(7.2) 



Im Folgenden werden die Gleichungen auf eine Dimension x beschrankt, da spater nur solche 
Systeme betrachtet werden. Fiir weitere Berechnungen ist es zusatzlich erforderlich, die Gr6i3e 
^ = in die x-Koordinate mit 

(7.3) 



X 

X ^ — 

£ 



aufzunehmen.^ 

Mit X andert sich Gl. ((7^ zu 



2m ^ 2m yHTg 
wobei hier zusatzlich das externe Potential 



m 



2m 

in die neuen Koordinaten eingesetzt wurde. 



-vl - Ngps. 



(7.4) 



(7.5) 



7.1.1 Die Phase der Storung 

Nehmen wir an, der Zustand sei 



ausniitzend, dass fiir eine Dimension 



Vs = —0x00 

m 



Ooix) 



m 



dx'vs{x') 



gilt. Die Phase von r} gibt die Realtivbewegung der Storung zum Kondensat an. 



(7.6) 



(7.7) 



^Die Terme ohne e konnen als Konstante zusammengefasst werden, sind damit fiir weitere Betrachtungen 
uninteressant. 

^Der Sinn dieser Tranformation zeigt sich spater, wenn die WKB-Naherung angewandt wurde. 
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Bevor dieser Ansatz in G1. (|7.4I) eingesetzt werden kann, wird er durch die neuen Koordinaten 

(Gi.dini)) 

ausgedriickt. 

Mit 'ip = £fi{t, x)e^~ ^ '"^^^ ^"^^ erhalt man somit eine Gleichung 

/ \ e^h^ ( \fp's'\ 

ihfi = Ngps i^fj + 7]^j - ie {2vs7]' + f/v'^) - - ~^^J ' C^-^) 

die nur noch vom Feldoperator fj abhangt. Setzt man in Gl. (|7.9|) eine Superposition zweier 
Feldoperatoren (0772 + 6%) ein, so muss jeder fiir sich diese erfiillen. Es liegt im Ortsraum 
keine Koppelung zwischen den Fluktuationen vor. Es bleibt zu untersuchen, wie die Situation 
im Impulsraum ist. 

Die Fourier-Entwicklung der Storung iiber einfache ebene-Wellen-Losungen 

Liegt der Feldoperator %l){x,t) vor, ermoglicht die diskrete Fourier- Analyse 

'iPk{t)=^ r'^'\-'''^u^,t) (7.10) 

und deren Riicktransformation 

vom Ortsraum in den Impulsraum zu wechseln und umgekehrt, wobei n G Z ist. 

Die Operatoren ipk{t) iiii A;-Raum konnen in Erzeugungs - und Vernichtungsoperatoren 
umgeschrieben werden: 

^k{t) = Uk{x)ak{t) + vl,{x)al^{t). (7.12) 

Im Heisenberg-Bild 

ak{t) = cike-'^^' (7.13) 
lasst sich der Feldoperator umschreiben in 



Mt) = E [M^)l^ke-'^'' + vl{x)ale'<'] . (7.14) 



45 



Mit dieser Entwicklung wird Gl. (|7.2() diagonalisiert, wobei die Eigenfrequenzen uj^ der Energie 

eines Quasiteilchens zugeordnet werden konnen. Dieses Verfahren wird als Bogoliubov- Transformation 

bezeichnet [H]. 

Die Amplituden und Wellenzahlen hangen hier vom Ort ab: 



erfiillen. Es ergibt sich, dass fiir ein Paar {uy^Vki) immer eine dazu redundante Losung 
(u^,M^) mit den Frequenzen ooy = — existiert. 

Die Entwicklung nach ebenen Wellen wird im nachsten Unterabschnitt in Gl. ()7.4() eingesetzt. 
Es ergibt sich jeweils eine Gleichung fiir die Erzeugung - bzw. Vernichtungsoperatoren, beide 
zusammen werden als Bogoliubov- Gleichung en bezeichnet. 

7.1.2 Die Bogoliubov-Gleichungen in einer Dimension 

Die Bogoliubov-Gleichung in ihrer bekanntesten Form ist 



vu(x)=v{x)e'^^^''-'^^-'^-'Trdx'v.(.') 



(7.15) 



Erzeuger und Vernichter miissen die Vertauschungsrelation 



(7.16) 




(7.17) 



mit 



ho{x) 



e^h^ 



+ Ve^t + 2Ngp, 



2m 



's 



und 



u,{x) = u{x)e'?^rdx'v.i.')^^rdx'k(.') 

v,{x) = ^{^^y-'7krd^'-s{.')^\r<^-'K-') 



Eingesetzt in die Bogoliubov-Gl. 1)7. 17|) . ergibt sich 



/i+n + NgpsV — ie 













(7.18) 
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wobei 

h± = —— + Ngps±h{kvo + ujj) . (7.19) 
ml 

eingefiihrt wurde. Mit diesen Gleichungen ist es moglich, die Dispersionsrelation uj{k) und 
die Amplituden der erzeugten Moden zu berechnen. 



7.2 Die Zeitentwicklung der Fluktuationen 

Ergeben sich aus den Bogoliubov-Gleichungen komplex imaginare oder negativ reelle Eigenfrequenzen, 
wird das Kondensat instabil. 

7.2.1 Dynamische Instabilitaten 

Die Zeitentwicklung ist durch Gl. (|7.13() gegeben. Fiir komplexe Frequenzen mit negativem 
Imaginarteil entwickeln sich die Moden exponentiell mit der Zeit {exp{uiit)). Der Erwartungswert, 
dass Moden an einem bestimmten Ort erzeugt werden, nimmt ebenfalls mit der Zeit exponentiell 
zu {exp{2u)it)) . Diese Zunahme erfolgt selir schnell. Die Fluktuationen konnen nicht mehr 
mit den linearisierten Gl. 1)7.2(1 beschrieben werden. 



7.2.2 Energetische Instabilitaten 

Fiir u; < befindet sich das System nicht im Grundzustand. In einem solchen System konnen 
eventuell Moden erzeugt werden, die Energie aus dem Kondensat befordern, wodurch der 
Zustand des Systems verandert wird. 
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Teil IV 

Schwarze Locher im 
Bose-Einstein-Kondensat 
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Kapitel 8 

Zigarrenformige Kondensate mit 
Senke 



Es wird ein eindimensionales Kondensat untersucht, dessen Bewegung auf die x-Achse eingeschrankt 
ist. Die Kondensatwolke nimmt eine zigarrenformige Form mit der Lange 2D an und wird 
zusatzlich symmetrisch zum Ursprung an zwei Stellen eingeschniirt. Dadurch ist im Bereich 
2L der Querschnitt der Kondensatwolke diinner als auBerhalb. Diese Konfiguration wurde 
im Artikel von Cirac et al. jl7j bereits behandelt. Die Aufgabe dieses Kapitels ist es, die dort 
prasentierten Ergebnisse eigenstandig zu berechnen. 



8.1 Das Modell 



In Abb. H8.1() ist das Dichteprofil der Kondensatwolke skizziert. Die Pfeile zeigen die Richtung 

der Hintergrundgeschwindigkeit an. Es sind zwei Strome eingezeichnet, die sich am Koordinatenursprung 

treffen. Es ist aus Griinden der Kontinuitat erforderlich, an dieser Stelle einen Laser anzubringen, 

der eine Senke bei {x = 0) erzeugt. Das Geschwindigkeitsprofil - der Gradient der Phase - 

soil innerhalb \x\ < L um einen Faktor a > 1 groBer sein, als auBerhalb \x\ > L + e. Ftir ein 

stationares Problem gilt mit G1. ()B.4|) . dass der Fluss pro Volumeneinheit und Zeit konstant 

sein muss, und damit 

c{x)'^v{x) = const. (8-1) 
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stationare 

Hintergrundgeschwindigkeit 




Abbildung 8.1: Schematische Darstellung des Dichteprofils eines zigarrenformigen Kondensats. 
Die Pfeile zeigen die stationare Hintergrundgeschwindigkeit des Kondensats v = h/m\/pQ an. Fiir 
Werte \x\ > L + e und \x\ < L ist die Fliefigcschwindigkeit konstant, wobei diese Konstante im 
letzteren Bereicli holier ist. Dieser Sachverhalt ist durch die Lange der Pfeile angedeutet. In den 
beiden Ubergangsbereichen L < \x\ < L + e verengt sich das Dichteprofil, so dass es, um die 
Kontinuitatsgleichung zu erfiillen, schneller fliefien muss. In der Mitte bei a; = trefTen die beiden 
Strome des Kondensats aufeinander, die von einem off-resonanten Laserstrahl aus dem Kondensat 
geleitet werden. 



Fiir die Schallgeschindigkeit wird 

c{x) 

gewahlt und liefert mit 



Co, \x\ < L 

co[l + {a - l)x/e], L<\x\<L + £ 
aco, L + e < \x\ 



v{x) 



die gewiinschte Verteilung fiir die stationare Hintergrundgeschwindigkeit. 



(8.2) 



.3) 



8.1.1 Die Gross-Pitaevskii-Gleichung fiir das zigarrenformige Kondensat 

Das Kondensat kann durch die zeitabhangige Gross-Pitaevskii-Gleichung 

(8.4) 
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d - - 
^^^^ = i^dx^^ + (A^o<7l^P + Ve^tix) - ivo6{x))^ 
at 2m 



beschrieben werden, wobei es erforderlich ist, den Zusatzterm —wqS^x) einzufiihren, der 
die Wirkung des Lasers an der Stelle x = beschreibt. Dort wird die Fliissigkeit mit der 
Geschwindigkeit vq aus dem Kondensat entfernt. 
Das Potential 

V,,t{x) = _ !!ly^ _ mc{xf. (8.5) 

ergibt sich durch Einsetzen des Grundzustandes <l>o in Gl. ()8.4|) . 



8.1.2 Stetigkeitsbedingungen am Punkt a: = 

Das Kondensat ist durch die Senke an der Stelle a; = in zwei Telle getrennt. Die Losungen 
miissen stetig ineinander iibergehen 

^(0+,t) - ^(0",t) = 0. (8.6) 

Es ist ^(0=^,t) = \\me-,Q^{±e,t). 

Eine Bedingung fiir die erste Ableitung ergibt sich durch Integration^ von Gl. (|8.4j) : 

^'(0+,t)-^''(0-,t) = -^^(0). (8.7) 



8.2 Die Wellengleichung 

Die Dynamik einer modulierten Dichtestorung wurde in Kapitel © hergeleitet. In diesem 
Abschnitt werden die Rechnungen auf das zigarrenformige Kondensat angewendet. Fiir ein 
symmetrisches System reicht es aus, das Verhalten fiir x > zu untersuchen. Am Punkt 
X = werden die rechten und linken Losungen stetig verbunden. 

^Hier wurde zusatzlich fiir beschrankte f{x) verwendet: 

lim / f{x)dx = ± lim f{Q)e = 

lim / 5{x')f{x')dx' = ±/(0) 
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8.2.1 Die Metrik fiir das zigarrenformige Kondensat 



Zur Berechnung der Wellenfunktion kann ^ <I) gesetzt werden, was bereits in Kapitel © 
berechnet wurde. Fiir eine Dimension wird Gl. (|l'{.12|) zu 

ihdMt,x) = (^-^dl + Vext{x) + NogMx)\^^ <^{t,x), (8.8) 

wobei der Term fiir die Senke nicht auftaucht, da nur Werte x > betrachtet werden. Das 
chemische Potential fi ist in V^xt enthalten. 

Induziert man von auBen eine lokale Dichteschwankung, breitet sich diese wie eine Schallwelle 
in einem Superfluid aus. Die Bewegungsgleichungen fiir die Storung ergeben sich, wenn in 
G1. H8.4|) fiir die Dichte p = po + epi und fiir die Phase 6 = 9q + eOi eingesetzt wird. Die 
Rechnungen entsprechen denen aus Kapitel Um Wiederholungen zu vermeiden, werden 
im Folgenden nur die Ergebnisse prasentiert. In erster Ordnung erhalt man fiir den Realteil 



Po = -—dxipodxOo), 
m 



i.9) 



und den Imaginarteil 



h 



Pi 



m 



dxipodxOi + PidxOo) 



Diese beiden Gleichungen lassen sich zusammenfassen zu einer Wellengleichung 

df.i^g'^'d.ei) = 0, 



(8.10) 



^.111 



mit der effektiven Metrik 
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-c^ + vl 
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-c2 
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^.12) 



Damit ist gezeigt, dass sich die Dichteschwankungen analog zu masselosen Teilchen 
in Anwesenheit eines auBeren Gravitationsfeldes ausbreiten. Als nachstes wird untersucht, 
welchem Gravitationsfeld die effektive Metrik entspricht. 
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8.2.2 Die kovariante Metrik 



Ein Umschreiben der Metrik (|8.12|) in ihre kovariante Form ergibt 







-Vx 








\ 




-Vx 
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/ 



.13) 



mit 500 = -c + v^. 

Ein schwarzes Loch liegt vor (siehe KapitelEj), wenn sich das Vorzeichen von qqq andert. 
Der Punkt mit g^Q = wird als Ereignishorizont bezeichnet. Damit kann iiberpriift werden, 
ob die efFektive Metrik mit der eines schwarzen Lochs vergleichbar ist. 



8.2.3 Der Ereignishorizont 

goo > 0: Nehmen wir an, fiir x < L sei > c{x)'^, dann gilt mit Gl. H8.2p und Gl.(jH3I), 

dass Vq > Cq, wobei vo = h/mdx0o{x) und cq = h/m \J A7rapo{x) sind. Da 
beide Geschwindigkeiten in diesen Bereichen konstant sind, kann man mit einer 
Konstante s > 1 schreiben, dass Vq = scq. 

goo < 0: Es muss moghch sein, die vorherige Bedingung und gleichzeitig v{x)'^ < c(x)^ fiir 
X > L + e zu erfiillen. In diesem Bereich ist v(x)'^ = v^/a^ und c(x)^ = c^Cg. 
Damit kann man erneut s = v{x)'^ /c{x)'^ = VQ/{cQa^) bestimmen. Fiir s < 1 ist 
dann f (x)^ < c(x)^ erfiiht, was durch geeignete Wahl von a erreicht werden kann. 

Es gibt einen Bereich, in dem die Hintergrundstromung schneller flieBt (supersonic) als 
die Storung sich ausbreiten kann. In diesem Bereich muss die Dichtemodulation in Richtung 
des Kondensats flieBen. Damit befindet sich diese im Inneren des schwarzen Lochs. 
Aufierhalb ist ein Propagieren in beide Richtungen moglich. Dort flieBt die Hintergrundstromung 
langsamer als die Schallgeschwindigkeit (subsonic). Dazwischen befindet sich der Ereignishorizont. 
Wenn e gegen Null geht, ist der Ereignishorizont bei xeh = L. 

^Es gilt, g^-'g^. = 1. 
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Geschwindigkeiten 




Ereignishorizont 



4 



3 



2 



1 



inner halb 



aufierhalb 







L 



X 



Abbildung 8.2: Fiir e ^ sind hicr die Quadrate der Schallgeschwindigkeit c und der lokalen 
Hintergrundgeschwindigkeit v in Abhangigkeit von x schematisch dargestellt. Es wurden hier cq = 1, 
s = 2 und (7 = 2 gewahlt. Innerhalb des schwarzen Lochs ist die Schallgeschwindigkeit kleiner, 
aufierhalb grofier als die Geschwindigkeit des Kondensats. 

8.3 Die Zeitentwicklung der Fluktuationen 

Aus Kapitel ist bekannt, dass Fluktuationen um den Grundzustand nach der Modenentwicklung 
die Bogoliubov-Gleichungen H7.18() erfiillen. Diese werden nun speziell auf das zigarrenformige 
Kondensat angewendet. 

Die Losungen der Bogoliubov-Gleichungen 

Es ist nicht notwendig, in den Bereichen < L und |x| > L + e die integrale Form ftir die 
Phase dx'v{x') zu verwenden, da dort die lokale Hintergrundgeschwindigkeit v konstant 
ist. Damit ist: 
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Damit lauten die Bogoliubov-Gleichungen (|7.18|) fiir den inneren und auBeren Bereich 
h^u + NgpsV — ie (^77— A; — h\v\^ u" + — ( — k"^ u =0 
h-V + NgpsU-ie \ —k + h\v\\ v" + '- [ —k" \ v =0 



.14) 



wobei 



f{2 ^2 

h± = + Ngps T {k\v\ + LOj) . 

ml 



(8.15) 



ist. Diese Gleichungen sind auf ganz M analytisch. Da c und v sich in den drei Bereichen - siehe 
Abb. H8.3() - unterschiedlich verhalten, werden fiir jeden Bereich die Bogoliubov-Gleichungen 
separat gelost und an den Grenzen verbunden. 



Die Randbedingungen 

Setzt man 'ip = efj{t, x)e'-^ ■f'' '"^'^">'^''' in die Stetigkeitsbedingungen Gl. (|Hli|) und Gl.lfHTjl ein, 
erhalt man fiir 77(0) 

7}(0+,t) -r}(0-,t) = (8.16) 

und 

7}'(0+,t)-r)'(0-,t) = 0, (8.17) 

nur gerade und ungerade Losungen (Abb.(jHiSI))- ist ausreichend, die rechte Seite x > 
zu untersuchen."^ 



Im offenen Intervall x > gehen an den Stellen x = L und x = L + e die Losungen der 
abgetrennten Bereiche ineinander liber. Wahrend innerhalb und aui3erhalb des schwarzen 
Lochs die Losungen explizit ermittelt werden konnen, kann durch die folgende Integralabschatzung 
fiir ein x£L<x<L + € 







J dx'v{x') 


<-i: 







dx 



[1 + ((7 - l)x/e]^ 



< voe < 1, 



der Zustand vereinfacht werden zu 



(8.18) 



Vu>,k 



Am Schlufi sind von alien Losungen nur die am Punkt x = Q stetigen erlaubt. 



.19) 
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wodurch die zweiten Ableitungen verschwinden. 



Das Verhalten an den Kontaktstellen wird von der Schallgeschwindigkeit bestimmt, deren 
zweite Ableitung 



c"(x) _ CT - 1 



6{\x\ -L)- -5{\x\ -L-e) 
a 



(8.20) 



c(x) e 

singulare Punkte hat.^ 

Analog zur Bestimmung von G1. H8.6|) und G1. H8.7|) lassen sich nun die Randbedingungen 
an den Stellen x = L und x = L + e herleiten: 

r)(L+)-7?(L-)=0 

r}'(L+)-f/(L-)=^r)(L) 

^ (8.21) 

r/(L + e+) -77(L + e^) =0 

rj'{L + e+) - fi'{L + e") =^^r/(L + e) 

ae 

In Abb. (|8.3j) sind alle Ergebnisse fiir die Randbedingungen zusammengefasst. Die beiden 
ebenen Wellen - die paarweise erzeugt werden - miissen jede fiir sich die Gl. (|8.21() erfiillen. 
Es ergeben sich zwei Gleichungen fiir Uk und Vk- 

Aus den Anschlussbedingungen und der Form der Moden fGl. (|8.19jl ) im Ubergangsbereich 
ergibt sich 

UkjL + e) _ VkjL + e) _ 
Uk{L) Vk{L) a' 

1st e klein, kann 

u{L + e+) ^ u{L+) + eu'{L+) 

als Taylor- Reihe bis zur ersten Ordnung entwickelt werden. Ein Umschreiben in u{L~^) = 
Uout{L) bzw. u{L~) = Uin{L) und Einsetzen in obige Gleichung ergibt zusammen mit den 



(8.22) 



'Zur Berechnung ist die erste Ableitung der Schallgeschwindigkeit Gl. 18.21 erforderlich: 

c'{L-)=0 c'{L+) = ^ 

c'((L + 6)-) = -^ c'((L + e) + )=0 
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Abbildung 8.3: Schematische Darstellung der Randbcdingungcn. Die schwarze Kurve zcigt eine 
Losung ftir .t > 0. Im Bereich L < x < L + e konncn die Wellenfunktioncn linearisiert werdcn und 
die drei Bereiche miissen stetig ineinander iibergehen. Eine Erweiterung der Losung auf die ganze 
x-Achse ist moglich. Die Stetigkeitsbedingung an der Stelle x = erlaubt nur gerade - griine Kurve 
- bzw. ungerade - rote Kurve - Wellenfunktionen. 



Anschlussb edingungen 



Uin,k = - eu'out kiL) + -Uout,k{L) 
a 

Vin,k = - e-Woui fc(^) + -Vout,k{^) 

a 



8.3.1 Zur Losung der Bogoliubov-Gleichungen 

Fiir Potentiale, die sich nur wenig mit der x-Koordinate andern, kann die sog. WKB- 

Naherung (siehe AnhangEJ verwendet werden. In diesem Fall sind die Schall- und Hintergrundgeschwindigkeit 

in der inneren und auBeren Region konstant. Mit Gl. (|7.5j) ist das externe Potential ebenfalls 

konstant. Zunachst werden mit der WKB-Naherung die Losungen von Gl. (|7.17() bestimmt. 

Anschliei3end werden mit den Randb edingungen die inneren und aui3eren Wellen miteinander 

verbunden. 
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Die WKB-Naherung 

Fiir ein sich wenig anderndes Potential lassen sich die Amplituden der ebenen Wellen als 
Taylor-Reihe schreiben: 

u{x) = uo + eui + .... (8.24) 
v{x) = vo + evi + .... (8.25) 

Eingesetzt in Gl. (|7.17|) ergibt sich in nuUter Ordnung in e das Gleichungssystem 

( ^+ ""^'Af ''']=0, (8.26) 

dessen Losungen die Determinante hj^h- — (mc^)'^ = verschwinden lassen. Das fiihrt zu 
einer Dispersionsrelation, 

+ (c2 - \v\^)k^ - 2\v\u;jk - = 0, (8.27) 

der en Losungen im Allgemeinen nicht zu finden sind, aber fiir spezielle lv analysiert werden 
konnen. 



Reelle Frequenzen Die Definition der Punktion 

/ Ji2 U2 

f^ = huJi = -hk\v\±hk\^— + 0"^, (8.28) 

V 4 

welche eine abgewandelte Form von G1. H8.27|) ist, erweist sich als sinnvoll. Die Interpretation 
dieser Formel wird in Abb. (|8.4j) erklart. Stellt man f± in Abhangigkeit von k dar - wie in 
Abb. (|8.4|) - ist es moglich, die Losungen von Gl. (|8.27j) graphisch zu ermitteln - vorausgesetzt, 
dass Lo reell ist. Die Losungen fiir eine spezielle Frequenz ergeben sich, wenn eine horizontale 
Gerade in die Abbildung eingezeichnet wird. Die Schnittpunkte mit den beiden Funktionen 
/+ und /_ entsprechen den Wellenzahlen der einzelnen Moden. In den verschiedenen Bereichen 
kann mit diesem Verfahren die Anzahl der reellen Wellenzahlen bestimmt werden. 

Zur Ermittlung der nicht rein reellen Wellenzahlen wird zusatzlich verwendet, dass Gl. (|8.27() 
ein komplexes Polynom vierten Grades ist und damit genau vier Nullstellen aufweist jl8j . 
Dieses Polynom besitzt nur reelle Koeffizienten, wodurch jede komplex Konjugierte einer 
Nullstelle ebenfalls eine Losung des Polynoms ist. Die Losungen treten in diesem Fall paarweise 
auf. 
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Abbildung 8.4: In den drei Plots ist die Energie f± - siehe G1. H8.28|1 - in Abhangigkeit von k 
aufgetragen. Bild (A) zeigt die Situation innerhalb des schwarzen Lochs. Dort treten fiir Energien, 
die unterhalb \E\ < huimax sind, vier reelle /c-Werte auf - oberhalb zwci. In (B) bzw. (C) - am 
Ereignishorizont bzw. aufierhalb des schwarzen Lochs - ist keine maximale Frequenz zu finden. 
Zusammenfassend ist zu vermerkcn, dass in alien drei Bereichen mindestens zwei Wellen reelle 
Wellenzahlen k mit cntgegcngesctzter Laufrichtung habcn. 



Innerhalb des schwarzen Lochs: < \v\^ 

Es existiert ein LVmax, bei welchem genau vier reelle Nullstellen in G1. H8.27() auftreten. 
Oberhalb dieser Prequenz liegen zwei, unterhalb vier reelle Nullstellen. Die Frequenz 
'^max isoliert die Bereiche verschiedener Losungen voneinander und wird deshalb als 
Bifurkationspunkt bezeichnet. 

Lo > uJmax ■ Zwei reelle und zwei komplexe kd und k^-^ Losungen. 

Lo = UJmax ■ Vier reelle Losungen, eine zweiter und zwei erster Ordnung. 

CO < Umax '■ Vier reelle Losungen. 

Am Ereignishorizont: = |fp 

Vw : Zwei reelle und zwei komplexe kd und k*-^^ Losungen. 

Aufierhalb des schwarzen Lochs: > 

Vu; : Zwei reelle und zwei komplexe kd und kd Losungen. 
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Rein imaginare Frequenzen: Einsetzen einer beliebigen, aber rein imaginaren Frequenz 
uj = ir in Gl.dini) liefert 



wobei fiir die Storung f]{x,t) = w{x)e eingesetzt wurde. 

Anstatt gleich nach den Moden zu entwickeln, wird zunachst die WKB-Naherung durchgefiihrt. 
In der nullten Ordnung in e erhalt man 



Diese Gleichung ist nur fiir wq = erfiillt. 

Zusammenfassend: uj kann nicht rein imaginar sein! 

Komplexe Frequenzen: Einsetzten einer beliebigen, aber festen komplexen Frequenz 
to = ujr + iLOi in Gl. ()8.27() liefert 



Nimmt man weiterhin an, dass k reell ist, dann konnen Real- und Imaginarteil dieser 
Gleichung voneinander separiert werden. 

Der Imaginarteil fiihrt zu einem Ausdruck fiir tUr = —kv, welcher zusammen mit dem reellen 
Anteil zu 



fiihrt, was fiir reelle k nicht losbar ist. 

Zusammenfassend: Fiir rein komplexe Frequenzen konnen die Losungen k nicht rein 
imaginar sein. 

Es bleibt zu klaren, wie die komplexen Losungen von k aussehen. Von ausschlaggebendem 
Interesse ist die Frage nach dem Vorzeichen des Imaginarteils. Mit Abschnitt (|8.3.1|) folgt, 
dass keine rein reellen Frequenzen k erlaubt sind. Das bedeutet aber auch eine Unterteilung 
der komplexen w-Ebene in zwei Bereiche. Der eine ist die positive komplexe Halbebene, 
der andere die negative. Der Grund hierfiir ist, dass alle Losungen zu einer Frequenz das 



ihTw{x) = Ngps{w{x) + w{x)'^) — i£{2vw{x)' + w{x)v') + 0{e^) 



ihVwo = Ngps{wQ + wl). 



(8.29) 




v^)k'^ — 2v{uJr + iuJi)k — uj^, + ujf — 2iuJrCVi = 0. 
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Vorzeichen ihres Imaginarteils ki nicht wechseln konnen. Ansonsten wiirde es beim Ubergang 
von +ki zu —kj einen Punkt geben, in dem ki = ware, was nicht erlaubt ist. Dadurch wird 
es moglich, kleine Prequenzen zu betrachten, deren Resultate dann auf alle to in der gleichen 
Halbebene erweitert werden konnen. 
G1. H8.27() reduziert sich ftir a; = zu 

k'l^^^ + {c'-v')^=0. (8.30) 
Als Losungen ergeben sich eine doppelte Nullstelle bei k = und zwei weitere Losungen mit 

Mit der Tatsache, dass fiir reelle uj < ujmax ahe vier Nullstellen ebenfalls reell sind, kann man 
k = ky. -\- ieki setzen, da der Imaginarteil verschwinden muss, wenn uji — > 0.^ 

Im nachsten Schritt werden die Annahmen - k = k^ + i£ki und uo = i£/[sec\ - in G1. H8.27|) 
eingesetzt 

^ {kt - e'^Qklkf + e^kf + i {e4kf.ki - eHkrkf)) + 
(c^ - v"^) (kl - e^kf + ie2krki) + e^2kiv/[sec] - ie2krV + e^/[secf = (8.31) 

und in Real- und Imaginarteil unterteilt. Weiterhin ist noch in die verschiedenen Ordnungen 
von e zu unterteilen. Es wurde beriicksichtigt, dass fiir u> ^ ie die Einheit der Frequenz 
l/[sec] nicht ohne Weiteres weggelassen werden darf. 



Reell und 0{e°) 



^,kt + {c^-v^)k^^=0 

kr = (8.32) 
k, = ±^y/v'^ -c^ (8.33) 



kr = in GLdHHU) 



Aw? 



e^kf + - (c2 - v^) £^kf/[sec] + £'^2kiv/[sec] + e^lsecf = 



Um Missverstandnisse zu vermeiden: Es ist moe; lich, in Gl.lESZJ 

rein imaginare Frequenzen einzusetzen, 
auch wenn diese in GL UT. 91 nicht erlaubt sind. Am Ende werden diese Losungen ausgeschlossen. Es ist nur 
ihr Verhalten in der oberen komplexen tj-Ebene von Interesse. 
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und 0(6^ 



-V ± c 



\sec\ 



(8.34) 



kr = ±^Vv2 - c2 in Gl.dHSU) 



h 



Am? " V h 

{v^ - c^) [-e'^k} ± ie^^fif^^k^^ + eHkil\sec\v =F ie^ x/u^TT^y + eVfsec]^ = 
und 0{e^) 



2 (z;2 _ 



sec 



In Tabelle 1)8. ip werden alle Ergebnisse zusammengefasst. 



k 

k — ~\~ 



ki -- 
k2 = 
fcs = + 
ki = — 



1 



sec 



sec 



C^ + ^T^Nc] 



+ i 



2 w2 

1 V 

2 t;2_c2 



sec 



Vorzeichen von /cj Vorzeichen von fcj 
fiir > fiir < c^ 



+ 
+ 



+ 



^.35) 



Tabelle 8.1: Fiir ein beliebig gewahltes komplexes to ist die Wellenzahl k ebenfalls komplex. Jeweils 
lassen sich zwei Losungen mit positivem bzw. negativem Imaginarteil ki finden. 



Zusammenfassend: Fiir kleine Frequenzen in der oberen Halbebene haben wir zwei 
Losungen mit positivem und zwei mit negativem Imaginarteil. Weiterhin haben Losungen 
ftir ujj > die gleiche Anzahl an positiven bzw. negativen kj. Das gilt fiir die obere w-Ebene 
und damit fiir die ganze Ebene. 

Da alle Ergebnisse auf der Giiltigkeit der WKB-Naherung beruhen, gilt es, diese zu iiberpriifen 

®Fiir v'^ < wird der Realteil von ks und ^4 imaginar, weshalb in diesem Fall die beiden zu verwenden 
sind. Daher sind fca = +i^\/v^ — und k4 = —i'^\J%P- — (? komlex-konjugiert zueinander. 
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Giiltigkeit der WKB-Naherung: Zur Kontrolle der WKB-Naherung ist zu analysieren, 
ob die Variation der Amplituden der ebenen Wellen klein ist. Dafiir ist es notwendig, uo{x) 
und vq{x) zu bestimmen. 

Dazu werden die Gl. H8.14p bis zur ersten Ordnung in e betrachtet. Es ergibt sich 



h+ui + NgpsVi - i 
h-Vi + NgpsUi — i 



—k{x) - hvj n'o + ^HxYuo 








^.36) 



Um die nullte Ordnung zu erfiillen, muss weiterhin die Gl. (|8.26() erfiillt werden, damit 
h^h^/Ngps = Ngps- Multiplizieren der ersten Zeile von G1. H8.36() mit h^/Ngps fiihrt zu 



NgpsUi + h^vi = i 



Ngps 



—k{x) -hv]uQ + —k{x)'uo 
2m J 2m 



Weiterfiihrend ist es notwendig, dass aus G1. H8.26|) 

uo h_ Ngps 



VQ 



Ngps 



''Const 



.37) 



folgt, was ^ = (loguo)' = (logt^o)' und h^/Ngps = hconst impliziert. Einsetzen in die zweite 



Zeile von Gl. lHM liefert 



const 



— -hv\ (log no)' + ^fc(^)' 



+ 



—k{x) + hv\ (logno)' + — M^)' 



0. 



Umformung und Integration kann gezeigt werden, dass uo{x) - und damit audi vq{x 
proprotional zu 



1 



duj/dk 

sind. Am Ereignishorizont wird der Nenner unendlich grofi, womit die Amplitude gegen Null 
gelit. Die WKB-Methode ist giiltig, da uq und vq am Ereignishorizont endlich sind. 



Auswertung der Randbedingungen fiir komplexe Frequenzen 

Zur Analyse dynamischer Instabilitaten sind nur die Losungen fiir komplexe Frequenzen 
von Interesse (siehe 1)7. 2() ). Es treten fiir ein beliebiges komplexes lo immer vier komplexe 
k auf, wobei zwei davon einen positiven und zwei einen negativen Imaginarteil haben. Fiir 
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Ini{k) < wachst die Storung exponentiell mit x (siehe Abschnitt 1)7. 2|) ) und entspricht 
einer raumlichen Instabilitat, wenn der Raum nicht begrenzt ist. Die auBeren Losungen sind 
daher auf Im{kout) > zu beschranken. 

Jede Linearkombination der ebenen Wellen in einem Bereich ist ebenfalls eine Losung. Alle 
Losungen miissen die Stetigkeitsbedingungen Gl. ()8.23|) erfiillen. 



An der Stelle x = L: Die inneren Losungen in allgemeiner Form sind 

—X^ TP .Jihn,j-Vo)(x-L) 

j 

,,. ('™^ _ \^ 77 u. .pi{kin,j+vo){x-L) 



a] 



wobei a = 1, 2 und j = 1, 2, 3, 4 ist J Gesucht ist die KoefRzientenmatrix 

/ Fii F21 \ 

F12 F22 

Fl3 F23 
y Fi4 F24 J 

und der Zusammenhang zwischen Uin und Vin, der aus G1. H8.37|) folgt: 

^aj — ^ ^m,l ^in,2 ^in,3 ^m,4 

Entsprechend ergibt sich fiir den auBeren Bereich 

„, ^ ('™^ _ pi(kout,m-vo)(x-L) 
"•out,a\-^) — / ^ ^ amy 

m 

„, ^ (^\ — \^ 77* i{kout,m+VQ){x-L) 

'-'out,a\'^) — / am.' ''QUI, j 
m 

wobei a = 1,2 und m = 1, 2 ist, und 

F* =( ^21 \ 

- U^, F^2 



und 



ho 



houtA hout,2 



(8.38) 



(8.39) 



.40) 



^.41) 



^.42) 



^.43) 



''Bei den KoefHzienten uq handelt es sich um Normierungsfaktoren, die als Eins gewahlt wurden. 
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Zusammen mit den Anschlussbedingungen (Gl.(jHI2Sl)) ergeben sich acht Gleichungen fiir 
zwolf Unbekannte Faj und F*^. Vier Unbekannte konnen daher frei gewahlt werden: 



1 
1 



.44) 



Mit dieser Wahl ergibt sich 



Uin,2 
Uout,l 
Uout,2 
Vin,l 
Vin,2 
Vout,l 
Vout,2 

wobei 



Fife. , (x— L) I 7-1 ifc. L) I j-i ik- „{x—L) , 7-1 ik .(x—L) 



_|_ _P22e ™.2(x-i,) _|_ _P23e 



+ F2Ae 



e 



/io«t,2e''^''-°-*-^("-^) 

(8.45) 



t2 



.46) 



ist. Einsetzen von Gl. ()8.45|) in G1. H8.23() und Umschreiben fiihrt zu 





1 


1 


1 


1 


\ 


/ 


Fal 


\ 






kin.2 


^m,3 


h~ 






Fa2 






hi 


h2 


00 








Fa3 




v 




h2ktn,2 


h3ktn,3 




/ 


V 


Fai 


/ 



( 



ok. 



out, a 



Auflosen nach Fq,j ergibt 



Fqj — M ^Ca, 



iA7) 
iA8) 
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mit 



M 



( 1 



1 



1 



m,l 



k 



in,2 

h2 



k 



in,3 

h3 



und 



( 



1 



(8.49) 



ak 



out.a 



V 



'i'(^hout,aki 

ah, 



out, a 

k: 



out, a 



outa,i^out,a 



i.50) 



I 



Es bleibt noch, die Losungen der rechten und linken Seite zu verbinden. 



Anschlussbedingungen an der Stelle x = 0: Die Anschlussbedingung Gl. (|8.6() an der 

Stelle X = lasst nur gerade f\g und ungerade r}„ zu. 



Ungerade Losungen: Fiir eine gerade Funktion gilt 
Einsetzen in G1. ()8.6I) liefert 



^.51) 



f?3(0,t) = -0. (8.52) 

Damit miissen Uin{0,t) und Vin{0,t) unabhangig voneinander an der Stelle x = 
verschwinden (tij„(0,t) = Vin{0,t) = 0). Da innerhalb zwei Losungen existieren, muss 
jede Linearkombination der beiden dort ebenfalls Null sein. Die beiden Gleichungen 

Uin,l{0) - Vin^l{0) = 
■Um,2(0) - Vin,2i0) = 



^.53) 



miissen gleichzeitig verschwinden. Zusammengefasst ergibt sich die Anschlussbedingung 
fiir gerade Losungen 



«m,l(0)-Um,2(0) - Uin,2iO)Vin,liO) = 0, 

in welche Gl. (|8.38|) eingesetzt werden kann: 



^.54) 



^.55) 
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Gerade Losungen: Fiir eine ungerade Punktion gilt 



f]u{x,t) = riu{x,t). 



(8.56) 



Einsetzen in Gl. (|8.6|l liefert 



7?„(0,t) + iVo7]u{0,t) = 0. 



(8.57) 



In Gl. (|8.38jl eingesetzt ergibt sich 




(8.58) 



8.3.2 Zur Losung der Gleichungen 

In der Veroffentlichung von Cirac et al. jl7j findet sich ein numerischer Losungsweg, der 
im Folgenden kurz skizziert wird. Die beiden Gleichungen H8.55() und H8.58() hangen fiir 
ein bestimmtes System nur von uj ab. Das System ist durch die Lange 2L der inneren 
Region, die Fliei3geschwindigkeit an der Senke vq, den Unterschied der inneren und auBeren 
Geschwindigkeiten a und durch die Konstante U = bestimmt. Es wird noch darauf 
hingewiesen, dass die Gesamtlange 2D der Kondensatwolke und die Grosse e des Ubergangsbereichs 
eine Rolle spielen, wobei hier nicht der exakte Wert ausschlaggebend ist. Ausreichend ist, 
das Kondensat moglichst groiB und die Verengung moglichst klein zu wahlen. 
Mit diesen Uberlegungen verbleiben zwei Gleichungen 



Zur Berechnung werden zunachst die SystemgroBen gewahlt und anschliefiend die Punktepaare 
{u, s) numerisch ermittelt, wobei beide im Allgemeinen komplexe Funktionen sind. Tragt man 
den Betrag von s iiber lo auf, reprasentieren die Nullstellen die Eigenfrequenzen. 

Interessant ist, mehrere Plots fiir unterschiedliche Langen des inneren Bereichs zu vergleichen. 
Es stellt sich heraus, dass die Anzahl der Instabilitaten empfindlich von L abhangt. Erst 
ab einem gewissen Duchmesser irh/kom — 5, wobei 5 sehr viel kleiner ist als vr/Zc, treten 
komplexe Eigenfrequenzen auf. Jede VergroBerung um tt/Zcq erhoht die Anzahl der komplexen 
Eigenfrequenzen um eins. Bei der auftretenden Wellenzahl handelt es sich um 



s{uj; a, U, vo,L) = 0. 




(8.59) 
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was physikalisch erklarbar ist: Erst wenn komplexe Eigenfrequenzen auftreten, werden Moden 
angeregt. Die Wellenzahlen dieser Moden sind in Tabelle (|8.1|) aufgelistet. Innerhalb des 
schwarzen Lochs kann es nur Losungen geben, wenn der Realteil ein Vielfaches von ^ ^ v^ — (? 
ist, da keine rein imaginaren fc-Werte erlaubt sind (siehe Kapitel H8.3.H) '). Die kleine Abweichung 
b entsteht, weil der Ereignishorizont keine exakt harte Wand ist. Damit konnen die dynamischen 
Instabilitaten auf die Bindungszustande zuriickgefiihrt werden. Erst wenn die Wellenlange 
der erzeugten Mode groB genug ist, um die Stetigkeitsbedingungen am Rand des Potentialtopfs 
zu erfiillen, konnen sich gebundene Zustande ausbilden. Damit gilt fiir die Frequenz des ersten 
gebundenen Zustands 

L Lko 

Obwohl mit der Lange 2L die Anzahl der gebundenen Zustande steigt, ergibt sich aus der 
Simulation fiir kleinere Locher eine hohere Instabilitat. 

8.4 Zusammenfassung 

Am Ende des Kapitels werden alle Ergebnisse zusammengefasst. Es wurden zwei Punkte 
ausgearbeitet. 

Zum einen betrachtete man eine von aufien angeregte Storung im Kondensat. Es resultierten 
dynamische Gleichungen fiir die Storung. Durch die Einfiihrung einer effektiven Metrik 
konnte gezeigt werden, dass sich die Phase der Storung wie Licht in der Nahe eines schwarzen 
Lochs bewegt. Die Eintrage in der Metrik bestehen nur aus Grossen des Kondensats. An dem 
Punkt, an dem die Hintergrundgeschwindigkeit des Kondensats dem Betrag nach gleich der 
Schallgeschwindigkeit ist, grenzen zwei Bereiche aneinander. In einem fliesst das Kondensat 
mit Uber schallgeschwindigkeit, im anderen unter Schallgeschwindigkeit. Im ersteren kann die 
Storung nur in eine Richtung propagieren, wahrend im anderen Bereich beide Bewegungsrichtungen 
moglich sind. Das zigarrenformige Kondensat wurde so gewahlt, dass dieses Verhalten im 
verengten Bereich auftritt. Gelangt die Storung in den eingeschrankten Bereich, wird sie 
vom Kondensat unaufhaltsam mitgerissen und am Koordinatenursprung iiber die Senke aus 
der Kondensatwolke ausgeschieden. Dieser Bereich entspricht dem Inneren eines schwarzen 
Lochs. 

Als Zweites wurde untersucht, wie sich der nichtkondensierte Anteil auf die makroskopische 
Wellenfunktion auswirkt. Fiir Temper aturen nahe am absoluten Nullpunkt der Temper atur 
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konnen die nichtkondensierten Atome als Storung im Bose-Einstein-Kondensat beschrieben 
werden. Diese regen Moden an, die als Quasiteilchen aufgefasst werden, es handelt sich um 
Phononen. Da nur dann Phononen produziert werden, wenn im Inneren des schwarzen Lochs 
gebundene Zustande auftreten, miissen diese dafiir verantwortlich sein. Je mehr gebundene 
Zustande moglich sind, desto mehr Moden mit komplexer Prequenz gibt es. Diese verursachen 
dynamische Instabihtaten, d. h. die Anzahl der Phononen wachst exponentiell mit der Zeit. 
Ein instabiles Kondensat verharrt nicht im Grundzustand. Die Rechnungen basierten auf 
der Annahme, dass es sich um kleine Storungen um den Grundzustand handelt. Damit sind 
die Bogoliubov-Gleichungen nur eine gewisse Zeit giiltig. 
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Kapitel 9 

Ringformige Kondensate mit 
Wirbel 

Als nachstes Modell wird ein Bose-Einstein-Kondensat betrachtet, in dessen Zentrum ein 
Wirbel ist. Es wird untersucht, ob in dieser Konfiguration ebenfalls ein schwarzes Loch 
auftritt und wie sich die Fluktuationen um den Grundzustand mit der Zeit entwickeln. Die 
Resultate dazu wurden bereits von Cirac et al. JJj verdffentlicht. Erneut ist es die Aufgabe 
gewesen, die Rechungen zu verifizieren. Neu in diesem Kapitel sind zwei Simulationen. In der 
Ersten wird gezeigt, wie sich eine ebene Welle im Kondensat verhalt. Das zweite Programm 
zeigt das Verhalten eines Wellenpakets in der Nahe eines Ereignishorizonts. 

9.1 Das Modell 

Ein Wirbel im Kondensat ist ein Punkt in der Kondensatwolke, um welchen das Kondensat 
rotiert. Die Rotationsgeschwindigkeit v (xl/r nimmt umgekehrt proportional zum Abstand 
zu. Die Dichte im Zentrum des Wirbels ist Null (siehe Abb. ((23)) ■ 

Die Bewegung ist in der r und z Achse durch ein externes Potential eingefroren (siehe 
Abb.dini)). 

Durch die Einschrankung der Bewegung auf eine Richtung kann der Zustand in Polarkoordinaten 

in 

$(t,r,7,z) =/(r,z)l>(t,7) (9.1) 
aufgeteilt werden. Beim Wechsel von kartesischen Koordinaten zu Polarkoordinaten andert 
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sich der Tangentialraum und damit der Gradient 





( 


dr 






r 1 




\ 





(9.2) 



und der Kotangentialraum. Damit ist bei der Integration zu beachten, dass 




dxdydz 




rdrd'fdz. 



(9.3) 



9.1.1 Die Gross-Pitaevskii-Gleichung fiir das ringformige Kondensat 

Wieder wird als Ausgangspunkt die zeitabhangige Gross-Pitaevskii-Gleichung Gl. (|3.12j) verwendet, 
die in den eingeschrankten Polarkoordinaten - dr = und dz = - umgeschrieben werden 
kann in 

ihdMt,r,-f,z) = (^-^^d^ + Ve.t{l) + Umt,r,j, z)\^^ ^t,r,j, z) , (9.4) 
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Abbildung 9.2: Es wird cin Kondcnsat (orange cingczcichnct) betrachtet, dcssen 
Hintergrundgcschwindigkeit v nur von 7 abhangt. Zur mathcmatischcn Bcschrcibung wcrdcn 
Polar koordinatcn vcrwcndet. 



wobei das chemische Potential bereits in Vext enthalten ist. Einsetzen von G1. H9.1|) und der 
dimensionslosen Zeit r = liefert 

idMt,^) = l-^d^ + VeM + ^mt,^)\A^{t,^), (9.5) 



mit dem dimensionslosen Potential 



Vext = VexJ^ (9.6) 



und der neuen Wechselwirkungskonstante 

U = U'^ J dzdrr\f{r,z)\'^, (9.7) 
wobei gleichzeitig der Grundzustand von $ auf die Gesamtteilchenzalil N 

(■2tx 

||.(7)|2 = Ar (9.8) 

^0 

normiert ist. 
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Der stationare Zustand fiir Gl. ()9.5|) ist: 



(9.9) 



9.1.2 Schall- und Hintergrundgeschwindigkeit im Kondensat 
Mit G1. (|4.H)) ist die Schallgeschwindigkeit im Kondensat durch 



c(7) = \/Up.,/N (9.10) 

gegeben und damit von der Dichte abhangig. Bei der Wahl der Dichte ist zu beriicksichtigen, 

dass im Bose-Einstein-Kondensat die Wirbel quantisiert sind. Integration der Hintergrundgeschwindigkeit 

entlang einer geschlossenen Bahn um den Wirbel ergibt 27rm, wobei m eine ganze Zahl ist 

und als Windungszahl bezeichnet wird. Die periodische Randbedingung im Kondensat kann 

mit der Windungszahl formuliert werden: 



1 



27r 



m = — d^vi^). (9.11) 

Als Wahl erweist sich 

N 

ph) = — {l + bcos{j)) (9.12) 

passend im Sinne der Randbedingungen und fiir die gewiinschte Analogie zum Schwarzen 
Loch, wobei < 6 < 1 ist. 

Zusammen mit Gl. (|4..Sj) und der Kontinuitatsgleichung fiir ein stationares Kondensat d.^{(?v) = 
sind die Schallgeschwindigkeit 



c{j) = J^{l + bcos{j)) (9.13) 



und die Hintergrundgeschwindigkeit 



l7rc(7j^ 



gegeben. 



9.2 Die Wellengleichung 

Die dynamischen Gleichungen fiir eine kleine Dichteschwankung im Kondensat - injiziert 
durch einen Laser - kann analog Kapitel (|8.2() hergeleitet werden. 
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9.2.1 Die Metrik fiir das ringformige Kondensat 



Ausgehend von der zeitabhangigen Gross-Pitaevskii-Gl. (|9.5|) wird fiir die Dichte p = po + epi 
und die Phase 6 = Oq + e6i eingesetzt. In erster Ordnung erhalt man fiir den Realteil 



drPi 



h 



m 



{pod^Oi + pid^Oo) 



und den Imaginarteil 



dr9i = -- {d^9oVei)-^pi] . 



(9.15) 



Diese beiden Gleichungen lassen sich in eine Gleichung fiir die Phase zusammenfassen 

- {ve[) - [vOi)' + ((c2 - v^) e[)' = 0, (9.16) 



die mit der effektiven Metrik 



-& + vt V. 



7 ^7 



\ 








1 
1 
1 



(9.17) 



zu 



-gg'^-'d^Oi) = 0, 



(9.18) 



umgeschrieben werden kann, mit v{t,0) = v{9). 

In Kapitel ()8.2() wurde anhand von qqq gezeigt, dass es sich bei Gl. ()9.17() um eine Metrik 
fiir ein Schwarzes Loch handelt. 



9.3 Die Nullgeodaten 

Die Nullgeodaten beschreiben die Bahnen der masselosen Teilchen im gekriimmten Raum. In 
diesem Fah ist das Linienelement Gl. (|5.6|) zeitartig {ds^ = 0). Speziell betrachten wir ebene 
Wellen 

^ = u;(7)e-^^*e^^'^T'=(^) (9.19) 
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die durch das Kondensat propagieren.^ Einsetzen in Gl. ()9.15() liefert 

= fw 

+ ie {w^' + 2w'^) (9.20) 

mit 

f = uj^ + 2vujk^ - {(? - v^)k^ (9.21) 

und 

1 = {(? - v'^)k - vuj. (9.22) 
Zur Losung dieser Gleichung wird die WKB-Methode verwendet (siehe Anhang Id))- 

9.3.1 Die WKB-Methode fur die Berechnung der Nullgeodaten 

Bei der WKB-Methode wird die Amplitude w{x) in eine Taylor-Reihe 

w{x) ^ wq + ewi (9.23) 

entwickelt, wobei fiir ein sich wenig anderndes Potential alle Terme mit hoherer Ordnung als 
0{e) vernachlassigt werden. 

Die Dispersionsrelation 

Fiir die nullte Ordnung muss die erste Zeile in verschwinden: 

f = Lo'^ + 2vujk^ - (c^ - t;2)A;2 = 0. (9.24) 
Fiir die Dispersionsrelation ergibt sich damit 

^(^) = ( f V (9-25) 
-v{-y) ± c(7) 



^ZuT Berechnung wurde 7 ^ im Exponenten eingefiihrt. Auf G1. H9.15|I angewendet ergibt: 
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Die Amplitude 

In erster Ordnung 

fwi + iiC'wo + 2^w'o) = 
ergibt sich die Amplitude wq, da / = gilt. Zu losen ist 

icwir = 

^^Wq = const. 



const. 



Mit Gl. (|9.25|) kann die Amplitude angegeben werden zu 

1 

Wo 



v{-i) - (-f (7) ± yc2(7)) 



Die Ebenen-Wellen im Kondensat haben die Form 



[t, 7) 



^ v{-i) - [-v{-i) ± \/c2(7) 



(9.26) 



(9.27) 



(9.28) 



9.3.2 Einlaufende und auslaufende Koordinaten 

Mit der Metrik G1. H9.17|) fiir das ringformige Kondensat kann das Linienelement 

ds'^ = -(c^ - vo^)dT'^ - 2vo^-i + 7^ , (9.29) 
(siehe G1. H5.6() 1 aufgestellt werden. Im Falle der Nullgeodaten ist ds"^ = und es gilt: 

= -(c-vo^,) - 2770^7 + 7^. (9.30) 
Die Losung der Gleichung lautet 

7 = —vq^ ± c (9.31) 
und durch deren Integration ergeben sich die einlaufenden 7* - und auslaufenden 7;!j_ Koordinaten 



7± 



/ 



-^(7') ± 0(7') 
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(9.32) 



Am Ereignishorizont v['^eh) = cilEH) divergiert j^. Ein Umschreiben des Linienelements 
in die neuen Koordinaten liefert 



ds^ = —c{(?duj^du-)^ 



(9.33) 



wobei zusatzlich u± = r — 7^ verwendet wurde. Je naher eine auslaufende Welle dem 
Ereignishorizont kommt, desto langsamer vergeht die Zeit. Es scheint fiir einen ruhenden 
Betrachter, welcher sicli aufierhalb des Bezugssystems der Welle befindet, als ob die Ausbreitungsgeschwindigkeit 
der Welle gegen Null geht, so dass diese das Schwarze Loch nicht verlassen kann. 

9.3.3 Simulation einer ebenen Welle im Kondensat 

Fiir u; = 10, 5 = 0.3, vc? = 1 sind in Abb. 1)9. 3() die Ergebnisse aus den vorherigen Abschnitten 
dargestellt. Die numerischen Berechnungen wurden mit MATLAB durchgefiihrt. Der Quellcode 
zum Programm befindet sich im Anhang (jD.lj) . 

9.3.4 Auswertung der Simulation 

In der Abb. ()9.3|l ist zu sehen, dass die ebene Welle, die sich mit dem Kondensat bewegt, 
keine Divergenzen hat. Die gegen den Strom des Fluids laufende Welle divergiert an zwei 
Stellen. Die Interpretation ist, dass beide Horizonte nur in gleicher Richtung zu iiberschreiten 
sind. Das bedeutet, dass hier zum schwarzen Loch noch eine weifie Quelle kommt. Zwischen 
den beiden (in Stromrichtung des Kondensats) befindet sich das schwarze Loch. Der erste 
Horizont (in Stromrichtung) kann nur in das schwarze Loch hinein passiert werden. Der 
zweite Horizont ist nur in Richtung aus dem schwarzen Loch hinaus zu liberqueren. Wahrend 
beim ersten das schwarze Loch nicht verlassen werden kann, ist es beim zweiten nicht moglich 
es zu betreten. 

9.3.5 Simulation eines Wellenpakets im Kondensat 

Mit einer Fourier- Analyse konnen die ebenen Wellen zu einem Wellenpaket iiberlagert werden. 
Damit ist es moglich, das Verhalten einer kleinen Dichtemodulation im Kondensat zu simulieren. 
Die Fourier- Transformation in diesem Fall ist 




(9.34) 
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Geschwindigkeiten 



2000 
1500 
1000 
500 


-500 
-1000 
-1500 
-2000 



Exponent f d-yk^ 



3.1416 

Winkel 7 



1.5708 3.1416 4.7124 

Winkel 7 

Einlaufende und auslaufende Ebene- Welle 





"0 

^^^^^ 
















1.5708 3.1416 4.7124 6.2832 

Winkel 7 - 
Ereignishorizont Ereignishorizont 

Abbildung 9.3: Im obcrstcn Bild sind die Schallgcschwindigkeit c und die stationare 
Hintergrundgeschwindigkeit v dargestellt. Im mittlcrcn Bcrcich flicfit das Kondcnsat mit 
Uberschallgeschwindigkeit, auBerhalb darunter. Es gibt damit zwci Schnittpunkte, die zwei 
Ereignishorizonten entsprechen. 

In der Mitte sind die Wellenzahlen aus G1. H9.28|I dargestellt. An den beiden Ereignishorizonten 
divergiert die Wellenzahl, wenn die Mode entgegengesetzt zum Kondensat propagiert. 
In der letzten Abbildung sind die einlaufenden und auslaufenden ebenen Wellen Gl. (|9.28() geplotet. 
Wahrend die auslaufende Welle an den beiden Ereignishorizonten divergiert, hat die Welle, die in 
das Schwarze Loch einlauft, kein divcrgcntes Verhalten in der Wellenzahl. Der Ausschnitt zeigt eine 
Vergroficrung der Wellen am Ereignishorizont. Es ist zu sehen, dass die auslaufende Welle - rot 
dargestellt - den Ereignishorizont nicht passieren kann. 
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Dichte 



T = tl T = 




1 7 
Ereignishorizont 



Abbildung 9.4: Zur Simulation wird cin Wellcnpakct bctrachtet, welches zum Zcitpunkt r = weit 
vom Ereignishorizont cntfcrnt ist. Dazu muss die Brcite des Wellenpakets klein sein im Vergleich zum 
Abstand des Ereignishorizont. 

wobei 

^^.(7) = «^o(7)e'^±" (9.35) 
ist. Die Fourier-Koeffizienten 0{u;) werden im Folgenden bestimmt. Zum Zeitpunkt t = wird 
ein GauB'sches Wellenpaket, weit entfernt vom Ereignishorizont betrachtet. Ist die Breite 
des Wellenpakets klein gegeniiber dem Abstand vom Ereignishorizont 70, konnen ^0(7) = 
fi^o(7o) und j^"-^ = 7^*'^'' (70) als Konstanten betrachtet werden (siehe Abb. ()9.4|l '). Die Fourier- 
KoefHzienten im flachen Raumen sind fiir ein GauB'sches Wellenpaket 

6(7) = e-^^^e^'^oT (9.36) 

Damit kann man die Fourier-Koeffizienten im Frequenzraum bestimmen: 

e{uJo) = Jki2^ [ 0(^)e-»'^o7i{7o)(7-7o)_ (9_37) 
w^o(7o)7r J 

Damit ist das Wellenpaket am Ort 70 zur Zeit r = bestimmt. Da die Fourier-Koeffizienten 
im Frequenz-Raum nicht vom Ort (bzw. Winkel) und der Zeit abhangen, bleiben diese fiir 
alle Zeiten und iiberall im Ort gleich. 

Zunachst kann man mit Gl. (|9.36() die Fourier-Koeffizienten G1. H9.37|) an der Stelle 70 be- 
rechen. Das Ergebnis wird dann in Gl. (|9.34|) eingesetzt. 

Die Simulation wurde wieder mit MATLAB durchgefiihrt, und der Quelltext befindet sich 
im Anhang (siehe (|D.2|) '). 
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Abbildung 9.5: Das erste Bild links oben zeigt die Welle bei r = 0, im zweiten Bild bei r = — 1 
und so weiter. Am Anfang handelt es sich um ein gaufiformig ausgeformtes Wellenpaket. Je langer 
es in der Zeit riickwarts propagiert, desto mehr verformt es sich. Den Ereignishorizont kann es nicht 
passieren. 



Auswertung der Simulation 

Das Ergebnis der Simulation ist in Abb. (|9. 5)1 dargestellt. Es ist zu sehen, dass es dem 
Wellenpaket nicht erlaubt ist, den Ereignishorizont riickwarts in der Zeit zu tlberschreiten. 
Anders formuliert, ein Wellenpaket, das sich auBerhalb eines schwarzen Lochs von diesem 
wegbewegt, kann nicht innerhalb des Ereignishorizonts gestartet sein. 

9.4 Die Strahlung 

Die Strahlung kann mit der Feldgleichung Gl. H7.4j) fiir das Kondensat berechnet werden, 
welche in Kapitel (|8.3|) bereits hergeleitet wurde. Fiir das ringformige Kondensat - mit der 
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dimensionslosen Zeit r und der auf eine Dimension eingeschrankten Bewegungsrichtung 7 - 
ist die Feldgleichung 



zl(r,7) = -i|(r,7)" + 



wobei U von Gl. (|9.7j) verwendet und zusatzlich 

vq = Rv 

eingefiihrt wurde. 

Der Zustands-Operator kann mit dem Bogoliubov-Ansatz entwickelt werden: 

Eingesetzt in Gl. (|9.38() ergibt sich eine Bestimmungsgleichung fiir fj: 

ir? = - ^ (^fl" - + i (^vov' + ^VQi]^ + ^po (rj + f?^) , 

9.4.1 Die Fourier- Transformierten Bogoliubov-Gleichungen 
In Kapitel (fTTTjl wurde gezeigt, dass es moglich ist, den Zustandsoperator fj 



(9.38) 



(9.39) 



(9.40) 



(9.41) 



nach ebenen Wellen zu entwickeln. Einsetzen in die G1. H9.4(J() fiir fj ergeben 



h—Uuj^n i 



(n - vo)u^,n - 2(^0^^^,") 



+ 



+ 



c 



+ CqV^^u =0 



wobei 



und 



h± = — + Cq± {nvo + u) 



U 



CO = ^PO 



(9.42) 



(9.43) 
(9.44) 
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eingefiihrt wurden. Ein Wechsel vom Orts- in den Impulsraum ermoglicht es, die beiden 
Gleichungen in der folgenden Form zu schreiben 



u, 



•Ju),n 



h+ f 

"■np Jnp 
~ fnp ^np 



"ui,p 



wobei 



1 

2^ 



fnp= — I d'ye 



2tt 



-i(n— p)7 



c(7)' 



und 



pMi) - TiMiY ^ { +-- 



lco(7)" 



2 CO (7) 



(9.45) 



(9.46) 



(9.47) 



sind.^ 

Diese beiden Fourier- Transformationen werden im Anhang exakt berechnet. Es ergibt 
sich: 



np 



2tt 



"n,p + —On,p+l + —On,p-l 



und 



htp = 7;{'^ + P) \/l - 62 as=p-n ± (fnp+ — -K.-D + 



4n2-U '^-b\ 

'n,p ~l~ „ Ps=n—p 



(9.48) 



(9.49) 



Die beiden Koeffizienten as=p-n und (3s=n-p sind im Anhang in G1. ()E.2() und G1. ()E.3|1 zu 
finden. Es handelt sich um zwei unendhche Reihen. 



9.4.2 Zur Losung der Gleichungen 

In der Veroffenthchung von Cirac et al. [21 hndet sich ein numerischer Losungsweg, der im 

Folgenden kurz skizziert wird. 

Bei groi3en Wellenlangen - kleinen n-Werten - treten Instabihtaten auf. Es kann gezeigt 

werden, dass die Einfliisse der hoheren Terme in den Fourier-KoefRzienten as=p-n und 

^Im Kapitel (18.311 wurde anstatt der Fourier- Transformation vom Orts- in den Impulsraum die WKB- 
Methode angewendet. Damit war es moglich, die zweiten Ableitungen von u und v zu vernachlassigen. Hier 
fallen die ersten und zweiten Ableitungen von und v^,ri durch die Transformation weg, da 

Mij,Ti(7)' = j d'yun^pe^''-""''^''^ = u„,p{-i{n - p))5„p ^ 
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(3g^n_p, auf die Instabilitaten vernachlassigbar sind. Es ist daher zulassig, einen Cutoff Q 
einzufiihren, der die Reilien der Fourier-Koeffizienten endlicli maclit. 

Die Bogoliubov-Gl.(jn321) enthalt nacli der Einfiilirung des Cutoffs 2{Q + 1) x 2{Q + 1) 
Matrizen, welclie numerisdi diagonalisiert werden konnen. Aus den sich ergebenden Moden 
sind jene erlaubt, welclie zusatzlich die Normierungsbedingung 



erfiillen. 

Die numerische Simulation zeigt, dass sowohl energetische wie dynamische Instabilitaten 
auftreten. 

9.5 Zusammenfassung 

Auch in diesem System ist es liber die effektive Metrik moglich, eine Analogie zum schwarzen 
Loch herzustellen. Eine von auf5en angeregte Storung hat die gleichen Bewegungsgleichungen, 
wie ein masseloses Teilchen im Gravitationsfeld. Speziell fiir lichtartige Teilchen, die als ebene 
Wellen beschrieben werden, liefert die Bewegungsgleichung die Nullgeodaten. Es ergeben sich 
zwei Losungen, die einlaufenden und auslaufenden Wellen. Erstere propagieren in Bewegungsrichtung 
des Fluids, die anderen dagegen. Die Simulation einer ebenen Welle mit einer festen Frequenz 
im Kondensat zeigt, dass an zwei Punkten Divergenzen fiir die auslaufenden Wellen auftreten. 
Zwischen diesen beiden Ereignishorizonten liegt (in Stromrichtung der Hintergrundgeschwindigkeit 
gesehen) das schwarze Loch. Innerhalb der Horizonte ist es nicht moglich, entgegen dem 
Fluidstrom zu propagieren. Der erste Horizont kann nur in Richtung des schwarzen Lochs 
iiberschritten werden, der zweite fiihrt unweigerlich aus diesem hinaus. Es ist nicht moglich, 
in das schwarze Loch iiber den zweiten Ereignishorizont zu gelangen, analog einer weissen 
Quelle im Universum. Uberlagert man alle Wellen zu einem Wellenpaket, kann das Verhalten 
der gesamten Storung simuliert werden. Es konnte gezeigt werden, dass ein Wellenpaket 
auBerhalb des schwarzen Lochs seinen Ursprung nicht innerhalb des Ereignishorizonts^ haben 
kann. 

Fiir die Geodaten ist vorausgesetzt worden, dass sich die Amplituden nur wenig andern. 
Die Amplituden der ebenen Wellen divergieren jedoch am Ereignishorizont. Im zweiten 
^Hier handelt es sich um den ersten Ereignishorizont, der nur Bewegungen in Richtung des schwarzen 




Lochs zulasst. 



86 







0.2 

b 



0.3 



0.4 



Abbildung 9.6: Darstcllung dcr dynamischen Instabilitatcn. In Abhangigkcit von den 
Kondcnsatgroficn {/ und h strahlt das System (weifie Bereiche), oder nicht (schwarze Bereiche). 
Kleinere Plots mit lioherer Auflosung zcigcn, dass cs sich bei den weifien Punkten tatsachlich um 
verbundene Bereiche handelt. Dem Artikel von Cirac et al. ^21 ist zusatzlich zu entnehmen, wie eine 
solche Konfiguration zu crzeugcn ist. 

Teil werden wieder die Fluktuationen um den Grundzustand behandelt. Hier treten die 
Divergenzen in den Amplituden nicht auf. Da diese Rechnungen um eine Ordnung genauer 
(im Bogoliubov-Ansatz) sind, haben die Divergenzen in den Amplituden keine physikalische 
Ursache. 

Die Bogoliubov-Gleichungen konnen durch eine Fourier-Analyse exakt berechnet werden. 
Als Losungsweg der sich ergebenden Gleichungen wurde von Cirac et al. ^2] ein numerisches 
Verfahren entwickelt. Zunachst wurden die Bogoliubov-Gleichungen vereinfacht und anschliefiend 
numerisch auf Diagonalform gebracht. Es zeigt sich in der Simulation, dass die Fluktuationen 
destabilisierend auf den Kondensatzustand wirken konnen. 

Es existieren stabile und instabile Bereiche. Ob ein System stabil ist oder nicht, hangt nur von 
den gewahlten Kondensatvariablen ab. Fiir die Instabilitat kommen sowohl energetische, als 
auch dynamische Ursachen in Frage. Die Zeitskala der energetischen [ixireeii < 0) Instabilitaten 
ist jedoch so grofi, dass sie vernachlassigbar sind. In Kapitel f\7.2\i ) wurde beschrieben, wie 
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die Fluktuationen sich energetisch instabil auf das Kondensat auswirken konnen. 

Nicht vernachlassigbar sind hingegen die dynamischen Instabilitaten, die fiir komplexe Frequenzen 

uj auftreten. Die numerischen Berechnungen ergeben, dass der Betrag der komplexen Eigenfrequenzen 

sensibel gegeniiber kleinen Anderungen der Kondensatgrofien (U, b, m*) ist. Fiir m* = 7 ist 

die maximale Anzahl der imaginaren Eigenfrequenzen in Abb. (|9.6|) dargestellt. 
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Teil V 

Das Universum im 
Bose-Einstein-Kondensat 
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Kapitel 10 

Das frei expandierende Kondensat 



Im Unterschied zu den beiden vorhergehenden Kapiteln wird in diesem Abschnitt ein Kondensat 
untersucht, das nicht nur vom Ort, sondern auch von der Zeit abhangt. Am Anfang ruht 
das Kondensat in einem aui3eren Potential, welches schlagartig abgeschaltet wird. Das nun 
freie Kondensat beginnt zu expandieren. Das Ziel dieses Abschnitts ist es, die Dynamik fiir 
eine von aufien injizierte Dichtemodulation zu berechnen. Es wird gezeigt, dass eine Analogie 
zwischen dem sich frei ausbreitenden Kondensat und der Expansion des de-Sitter-Universums 
besteht.^ 

10.1 Das Modell 

An das aufiere Potential wird nur die Bedingung gekniipft, dass es spharisch und dreidimensional 
sein muss. Dadurch wird erreicht, dass nach dem Abschalten des Potentials die Kondensatwolke 
homogen radial expandiert. Punkte im Kondensat bewegen sich auf radialen Bahnen, so dass 

r{t) = z{r,t = 0)b{t) (10.1) 

gilt (siehe Abb. (|lU.l() ). Die neue Koordinate z ist zeitunabhangig. Mit der Produktregel fiir 
die Ableitung von r{z,b) ergibt sich die radiale Hintergrundgeschwindigkeit 

Vr = Ir. (10.2) 



^Erwahnenswert ist eine bisher unverolTentlichte Arbeit von L. Garay 1101 . welche sich ebenfalls damit 
beschaftigt. 
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Kondensat zum Zeitpunkt t = 
Kondcnsat zum Zeitpunkt t — ti 




Abbildung 10.1: Zu Zeiten t < befindet sich das Kondensat (orange) in einem spharischen 
Potential. Am Zeitpunkt t = verschwindet das Potential durch pl5tzliches Abschaltcn. Die 
Kondensatwolke beginnt spharisch zu expandieren. 



Das Volumen der Kondensatwolke am Zeitpunkt t = ist 

Vo = /<fa<i,d. .(,-,* = 0). 
Mit der Zeit wachst das Volumen auf 

V{t) = b{tfVo 

an. Das Volumen ist indirekt proportional zur Dichte. Fiir Gl. (|4.3() ergibt sich 

1 



c{t,r) 



co{r), 



(10.3) 



die radiale Schallgeschwindigkeit. Ziel ist es, die Funktionen fiir c{t, r) und v{t, r) derart zu 
bestimmen, dass eine Analogic zum de-Sitter-Universum (siehe Absclinitt (|5.1.2|) ) besteht. 



10.2 Die Wellengleichung 

Im Falle eines frei expandierenden Kondensats ist Vext{t > 0,r) = 0. 

Induziert man mit dem Laser zusatzlich eine Storung in der Kondensatwolke, werden die 
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Dichte 



und Phase 



p{f, t) = po{f, t) + £pi{r, t) 
9{f,t) = 9o{f,t)+eeiif,t), 



minimal moduliert. 

Die Berechnungen fiir die Metrik sind in den Kapiteln © und zu finden. Es handelt sich 
hierbei nicht um exakte Kalkulationen. Es war erforderlich, die Thomas-Fermi-Naherung 
(G1. ()B.1|) ') zu verwenden. Diese gilt aber nur, wenn die kinetische Energie klein im Vergleicli 
zur restlichen Energie ist: 

Exin Epot + Eint- 

Am Anfang der Expansion ist 

ERin ^ Epot 

diese Bedingung gewahrleistet. Mit der Zeit wandelt sich potentielle in kinetsiche Energie 
um. Es ist notwendig, die Rechnungen zeitlich einzuschranken. Fiir die Storung muss - um 
die Thomas-Fermi-Naherung anfangs verwenden zu konnen - zusatzlich gelten: 



< Ngpi 
AmNg 



(10.4) 



Unter diesen Voraussetzungen ergibt sich fiir die Metrik: 
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Damit kann die Wellengleichung geschrieben werden als 

was einer Wellengleichung mit Metrik und damit im gekriimmten Raum entspricht. 

Der Unterschied zu den beiden vorherigen Kapiteln ist, dass es sich hier um Geschwindigkeiten 
v{t, r) und c(t, r) handelt, die explizit von der Zeit abhangen. Die Herleitung in Kapitel © 
ist jedoch allgemein gehalten und nicht auf ein zeitunabhangiges Kondensat beschrankt. 
Als nachstes wird die Metrik an das frei expandierende Kondensat angepasst. 
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10.3 Das Linienelement fiir ein spharisches Kondensat 



Die Kovariante zu g^^ ist 



\ 



c2 - V^) 


-Vx 


-Vy 


— V 


-Vx 


1 








-Vy 





1 





-Vz 








1 



Das Linienelement aus G1. (|5.6I) fiir diese Metrik ist 

df = c [-(c^ - v'^)df - 2vdfdt + dx^ + dy^ + dz^] 

wobei 

dr = 



(10.5) 




verwendet wurde. Da die Kondensatwolke radial expandiert, ist es sinnvoU, von den bisher 
verwendeten kartesischen in spharische Koordinaten zu wechseln. 

Fiir eine radiale Hintergrund - und Schallgeschwindigkeit kann das Linienelement in Kugelkoordinaten 
angegeben werden: 



^ [-{c{r,tf -v{r,tf)dt'^ - 2v{r,t) drdt + dr"^ + d^"^] 



(10.6) 



Um das Linienelement zu diagonalisieren, wird r in den neuen Koordinaten ausgedriickt. 
Dazu wird das erste 



dr = b{t)dz + ^r{t)dt 

0[t) 



(10.7) 



und zweite 



dr^ = h^dz^ + 2bvrdzdt + v'^dt'^ 
totale Differential berechnet. In Gl. (|10.5)) eingesetzt ergibt sich 

= coir) \-co{t)^ b{ty'/^ dt" + b{tfl^ dz" + 6(t)V2 ^2 ^^2 



(10.. 
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10.4 Die de-Sitter-Metrik im Kondensat 



In der de-Sitter-Metrik hangt goo weder von der Zeit, noch vom Ort ab (siehe G1. H5.1U|) ). Die 
Zeitabhangigkeit kann durch die Transformation 



7" = y ^/goodt 
dr =h{t)-^/^dt 



(10.9) 



beseitigt werden. Die Ortsabhangigkeit kann nicht beseitigt werden. 

Fiir ein konkretes auBeres Potential kann gepriift werden, ob ein Bereich existiert, in 
dem die Dichte annahernd konstant ist. Zum Beispiel ist der Grundzustand im harmonischen 
Potential zum Zeitpunkt t = durch 



$o(t = 0,f) 



'/i — i'mu!{t)^f^ 



Ng 



;io.io) 



gegeben. Hierbei handelt es sich um die Kugelgleichung. In Abb. (|10.H) befand sich das 
Kondensat in einem dreidimensionalen harmonischen Potential. Eine Einschrankung der 
Beobachtung auf den Bereich um das Zentrum des Potentials erlaubt es, die Schallgeschwindigkeit 
als konstant zu betrachten. 

Das Linienelement mit der neuen Zeitkoordinate r auf einem beschrankten Bereich im 
Kondensat ist damit: 



= coir) l-coirfdr^ + b (t(r))^/2 dz^ + b (t(r))^/2 z'^dn'^ 



(10.11) 



Ein Vergleich mit der de-Sitter-Metrik zeigt, dass 



/6(t) = exp(2iJr) 



(10.12) 



^Der Grundzustand kann berechnet werden, indem in Gl. 13.811 das harmonische Potential Vext = ^mLo{t)^f'^ 
eingesetzt wird. Mit der Thomas-Fermi-Naherung kann die kinetische Energie vernachlassigt werden. 
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sein muss (siehe Abschnitt 1)5. 10() ). Durch Integration von G1. ()1U.9|) 

r dr' = f h{t')dt' 
Jo Jo 

r dT' = [ exp{AHT)-^/Ut' 
Jo Jo 

[ dt' = I expiHrfdr' 
Jo Jo 

t = — exp(9ifr) — 

^exp(9ifT) = 9Ht + l 
kann 6(t) in Abhangigkeit von t angegeben werden (siehe jlUj'): 

6(t) = (9m + 1)^/9 (10.13) 

In der VerofFentlichung von Castin et al. [20] wird mit einer Skalierungs- Transformation 
die Zeitentwicklung fiir ein frei expandierendes Kondensat berechnet. Das Kondensat befindet 
sich vor der Expansion in einem dreidimensionalen harmonischen Potential. Die Formel fiir 
Lo in bewegten Koordinaten ist 

m = ^,-io{tfb{t). (10.14) 



In |lUj ist die Frequenz ebenfalls berechnet worden. Es ergibt sich: 

coit) = uoHt)-'/' (l + ^6(t)i/2) . (10.15) 

10.5 Zusammenfassung 

Fiir das frei expandierende Kondensat kann eine effektive de-Sitter- Metrik eingefiihrt werden. 
Das Problem der ortsabhangigen Schallgeschwindigkeit lafit sich durch Einschrankung auf 
einen kleinen Raumbereich, in dem die Dichte hinreichend konstant ist, beseitigen. Befindet 
sich das Kondensat vor Beginn der Expansion im harmonischen Potential, werden nur Storungen 
im Zentrum des Potentials betrachtet. Wahrend der Ausbreitung nimmt die kinetische 
Energie zu. Fiir kleine Storungen (siehe G1. H10.4|) ) ist die Thomas- Fermi-Naherung dennoch 
giiltig, solange die Energie der Storung klein gegeniiber der Wechselwirkungsenergie ist. Mit 
der Grofie der Kondensatwolke nehmen die atomaren Abstande im Kondensat zu, wodurch 
die Wechselwirkung abnimmt. Die Gleichungen werden ungiiltig. 
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Teil VI 

Zusammenfassung und Ausblick 
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In alien drei Kondensaten - zigarrenformig, ringformig und frei expandierend - konnte eine 
effektive Metrik eingefiihrt werden. Das zigarren - und ringformige Kondensat sind geeignete 
Systeme zur Simulation eines Schwarzen Lochs. Eine injizierte Dichtemodulation, die gegen 
die Stromrichtung des Fluids propagiert, divergiert am Ereignishorizont. Dieser Effekt kann 
im Experiment nachgewiesen werden. Ohne auBere Einwirkung werden Phononen emittiert. 
Diese Strahlung entsteht aus dynamischen Instabilitaten, die durch gebundene Zustande 
im Inneren des Schwarzen Lochs hervorgerufen werden. Die hier verwendete Bogoliubov- 
Naherung verliert jedoch nach kurzer Zeit ihre Giiltigkeit. Das Zeitfenster fiir die Messung 
der Strahlung ist daher nicht unbegrenzt. Die Beantwortung der Frage, ob es sich bei dieser 
Strahlung um Hawking-Strahlung handelt, ist schwierig. Zum einen bleibt offen, ob es sich 
um lokale, oder globale Strahlung handelt. Die Hawking-Strahlung eines Schwarzen Lochs 
ist lokal an den Ereignishorizont gebunden. Die angeregten Moden im Kondensat treten 
liber all auf. Es bleibt aber zu priifen, ob die Amplitude der Strahlung vom Ort abhangt. Sie 
miisste iiberall, auBer am Ereignishorizont, verschwindend klein sein. Komplizierter ist das 
Problem der Stabilitat. Schwarze Locher werden in der Gravitationsphysik als weitgehend 
stabil betrachtet. Instabil erzeugte Teilchen konnen als kleine Korrektur zur eigentlichen 
Strahlung betrachtet werden, d. h. die Zeitskala ist groB gegeniiber der Zeitskala fiir stabile 
Hawking-Strahlung. 

Das frei expandierende Kondensat kann als dynamisches Modell eines flachen Universums 
verwendet werden. Es wurde die Metrik fiir eine sich radial ausbreitende Kondensatwolke 
bestimmt. Unter der Bedingung, dass die Dichte in einem hinreichend grofien Bereich annahernd 
konstant ist, und einer Einschrankung in der Zeitdauer, ist es theoretisch moglich das de- 
Sitter-Universum im Bose-Einstein-Kondensat zu simulieren. In der Arbeit von Garay ^0] 
wird gezeigt, dass die Analogic zum de-Sitter-Universum erst ab drei Dimensionen moglich 
ist. Die Frage, ob im frei expandierenden Kondensat Phononen erzeugt werden, bleibt offen. 
In diesem Zusammenhang konnte versucht werden, dem Spektrum der Quasiteilchen eine 
Temperatur zuzuweisen. Mit einem solchen Modell konnten verschiedene Szenarien fiir unser 
Universum im Labor getestet werden. Ahnliche Anstrengungen wurden in |^ verfolgt. Dort 
wurde fiir ein oszillierendes Kondensat die Strahlung berechnet. Es stellte sich her aus, dass 
ein zweidimensionales Kondensat stabil ist, wahrend im dreidimensionalen Fall Phononen 
erzeugt werden. 
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Teil VII 



Anhang 
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Anhang A 

Quantenmechanische Streuprozesse 



Quantenmechanische Streuprozesse unterscheiden sich von klassischen Streuprozessen insofern, 
dass der Prozess nicht deterministisch ist. Es konnen nur Wahrscheinlichkeiten fiir die einzelnen 
Streuwinkel angegeben werden. 

A.l Reduziertes Zwei- Atom-Problem 

Sind zwei Atome durch ein Potential V{fi2) verbunden, wobei ri2 der Abstand ist, kann das 
System zu einem Ein-Teilchen-Problem reduziert werden. Die Bewegung im Schwerpunktsystem 



Abbildung A.l: Die Bewegung zweier wechselwirkender Atome zueinander kann auf die Bewegung 
eines Tcilchcns um den Schwerpunkt des Systems mit der reduzierten Masse fi — mim2/mi +TO2, 
im aui3eren Potential beschrieben werden. 




zwei Teilchen 

zwei interne Wechselwirkungspotentiale 



ein effektives Teilchen 
ein aufieres Potential 
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wird bestimmt durch 

Hrel = ^Vl,+Vin2), (A.l) 

wobei hier mi = m2 = m, und damit die Masse des Schwerpunkts /i = m/2, fi2 = r2 — ^i, 
der Verbindungsvektor und der Relativimpuls pi2 = P2 — Pi als Ausgangspunkt ftir das 
Korrespondenzprinzip^ verwendet wurde. 



A. 2 Streuung zweier Atome 

Um die Streuung zweier Atome zu beschreiben, bleiben wir in dem zuvor eingefiihrten 
Schwerpunktsystem, in welchem die Relativbewegung durch die Dynamik eines Zustandes 
beschrieben werden kann. Fiir ein zeitunabhangiges Potential sei der gestreute stationare 
Zustand ip{r'i2) ein Eigenzustand zum relativen Hamilton-Operator Hrei mit den positiven 
Eigenwerten E = h'^k^/2jj,. Die Schrodinger-Gleichung fiir das Problem lasst sich dann 
schreiben als 

(A + k^) ^(fi2) = '^Vifumfu). (A.2) 

Zusatzlich zu den gestreuten Zustanden konnen auch gebundene Zustande auftreten. Zur 
Losung dieser inhomogenen DifFerentialgleichung ist die integrale Form von G1. (|A.2|) 

(A + k^) i;{n2) = J d'ri2^V{f'i2)Hri2)^in2 - f^)- (A.3) 

aufschlussreich. Dies ist eine Differentialgleichung mit einer ^-Distribution als Inhomogenitat. 
Die Losungen solcher Differentialgleichungen sind bekannt, sie konnen mit Hilfe von sogenannten 
Greenschen Funktionen direkt angegeben werden. Die 5 — Distribution auf der rechten Seite 
von G1. ()A.3|) wird dazu durch 

(A + k^) i^aifu) = S{n2 - f'u) (A.4) 
ersetzt. Man bezeichnet die Losung tpQ = — ^^-pr^— als Greensche Funktion zum Differentialoperator 

{A + e). 

Um die vollstandige Losung ^(ri2) zu erhalten, muss zuvor noch die Losung der homogenen 
Differentialgleichung ipo{ri2) addiert werden. Fiir diese muss 

(A + k^) Mru) = (A.5) 

^Das Korrespondenzprinzip ersetzt in der Hamilton- Funktion Observablen durch Operatoren, hier Impuls 
durch Impulsoperator pi2 — > fVp^j. 
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gelten. Der gestreute Zustand ist dann 

W12) = Mri2) - ^ / ^M2^?^^^("i2)V'(rl2)- (A.6) 

Bei der homogenen Losung "00(^12) = exp{ikr 12) handelt es sich um eine einlaufende freie 
Welle, wobei die partikulare Losung eine auslaufende Kugelwelle darstellt. 
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Anhang B 



Berechnungen der 
hydrodynamischen Gleichungen 

B.l Berechnung der einzelnen Terme 

Hier werden die Terme aus Gl. ()6.1() 



+ ( Nog I Vpo + epi I ) V Po + epie 



berechnet. 



ihdt^/poTTpie 



i(6»o+£0i) 



„i(6»o+e6li) 



gi(6»o+£:0i) 



-ehOi 

h 
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„i(6»o+e6ii) 



gi(6»o+e6»i) 
+i , =X 



+e— V^oV^i 
m 

|;V(«V«„) 



-e—V{piV9o + poVei) 
2m 



In den Rechnungen konnte der Term 



vV(/oo + gpi; 



mit der Thomas- Fermi- Ndherung 



max Koo 



(B.l) 



vernachlassigt werden. 
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+eNQgpi] 

B.2 Die Kontinuitats- und Hamilton-Jakobi-Gleichung. 

Die Gleichungen werden einerseits nach Real - und Imaginarteil und andererseits in die 
verschiedenen Ordnungen aufgeteilt. 



Realteil mit 

e, = -^(vo,f-J--^ (B.2) 

2m Vext n 

0i = --{V0oy0i)-^Pi (B.3) 
m a 

Po = --V(poV^o) (B.4) 
m 

Pi = --V{poVei + piV6o) (B.5) 
m 

Gl. HB.2|1 und Gl. HB.3|) entsprechen dem quantenmechanischen Pendant der Hamilton-Jakobi- 
Gleichung in der klassischen Mechanik. Die beiden letzten Gleichungen sind die Kontinuitatsgleichungen 
fiir das Kondensat G1. HB.4() und G1. HB.5() fiir die Storung. 



Realteil mit 



Imaginarteil mit 



Imaginarteil mit 
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Anhang C 



WKB (Wentzel-Kramers- 
Brillouin)-Methode 



Die ausfiihrliche Behandlung der WKB-Methode wird motiviert durch die haufige Anwendung 
in dieser Arbeit. Der Grundgedanke ist, dass ein sich wenig anderndes Potential den Erwartungswert 
- eines sich im Potential bewegenden Teilchens - wenig beeinflusst. Die Amplitude der 
Wellenfunktion kann in einer Taylor- Reihe entwickelt werden, wobei dessen Glieder in hoherer 
Ordnung verschwinden. Dieser Zusammenhang wird im Folgenden ausfiihrlich erortert. 

C.l Aufenthaltswahrscheinlichkeit in einem auBeren Potential 

Propagiert eine Welle in einem aufieren Potential, andert sich die Wellenfunktion. 

C.1.1 Klassische Wellenfunktion in einem aui3eren Potential 

Die klassische Wellenfunktion wird beschrieben durch 



ct^ui{x)=we-'^'e- 



(C.l) 



Im konstanten Potential ist die Aufenthaltswahrscheinlichkeit 




(C.2) 



iiberall gleich grofi fAbb. (|C.l|) ). 



Ill 




x 



V{x) 



X 




'X 



Abbildung C.l: Die linke Halfte stellt cine Welle (schwarz) in einem konstantcn auficrcn Potential 
(griin) dar. Die Aufenthaltswahrscheinlichkeit, das Teilchen in einem bestimmten Intervall zu finden, 
hangt nicht von dessen Lage, sondern nur von dessen Lange ab. Fiir ein ortsabhangiges Potential 
(rechte Seite) ist die Aufenthaltswahrscheinlichkeit nicht mehr unabhangig von x. Nahe an einem 
Minimum ist es wahrscheinlicher das Teilchen zu finden, als bei einem Maximum. 



Im Potentialtopf wird der Erwartungswert - Amplitude w = w{x) - ortsabhangig, wahrend 
die Wellenzahl k unverandert bleibt. 



C.l. 2 Quantenmechanischer Zustand in einem aufieren Potential 

Die quantenmechanische Beschreibung eines Teilchens 



gleicht der ebenen Welle, jedoch mit dem Unterschied, dass zusatzlich die Wellenzahl k = 
k{x) vom Ort abhangt, wenn es sich in einem auBeren Potential befindet fAbb. (|C.2|) ). 




(C.3) 
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V [x] — 




X 



Abbildung C.2: Fiir cin variicrcndcs Potential (griin) andcrt sich auch dcr Zustand (schwarz). Die 
Wcllcnzahl und die Amplitude sind ortsabhangig. 



C.2 Die WKB-Naherung 

1st die Anderung mit dem Ort schwach, kann in G1. ()C.3|) die Amplitude w{x) durch 



w{x) = wq + ewi + ... - 



(C.4) 



eine Taylor-Reihe entwickelt werden, vorausgesetzt, die Terme hoherer Ordnung konnen 
vernachlassigt werden. 



C.3 Giiltigkeit der WKB-Naherung 

Die WKB-Naherung ist giiltig, wenn wi <C wq ist. In der Talyor-Reihe sind die Terme der 
Ordnung n proportional zur nten Ableitung. Die erste Ableitung der Amplitude nach dem 
Ort ist 

d,^{x) = + ikix)) tl;{x), (C.5) 
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damit fiir ^ 

<.k, (C.6) 



w 

d. h., wenn sich die Amplitude nur wenig mit dem Ort andert. 
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Anhang D 



Quelltext fiir Programme 

D.l Programm zur Simulation einer Ebenen- Welle im ringformigen 
Kondensat 



°/o Festlegung der Konstanten 



anzahl=10000; 
const=l ; 
xanf aiig=0 ; 
xende=2*pi ; 
omega=10; 
b=0.3; 



% Anzahl der diskreten Schritte. 
7o Es gilt: v_s*c~2=const 



7, Frequenz 

°/o Konstante fiir die Dichtefunktion. Es gilt: 0<b<l 



% Hauptteil 

x=linspace (xanf ang,xende, anzahl) ; 
dichte_s=(l+b*cos(x)) ; 
v_s= (const . /dichte_s) ; 



% x-Werte 
% Dichte 

% Hintergrundgeschwindigkeit 
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7o Berechnung der Amplituden: 

f aktor (1 , : )=1 . /sqrt (abs (v_s-(-v_s+sqrt (dichte_s) ) *omega) ) ; % fur entgegenlauf ende Welle 
faktor(2, : )=1 ./sqrt(abs(v_s-(-v_s-sqrt(dichte_s))*omega)) ; % fiir mitlaufende Welle 

% Berechnung der Koordinaten: 

w(l,:) =1 . /(-v_s+sqrt (dichte_s) ) ; % fiir entgegenlauf ende Welle 

w(2,:) =1 . / (-1) * (-v_s-sqrt (dichte_s) ) ; % fiir mitlaufende Welle 

7o Berechnung der Ereignishorizonte : 
EHl=acos( (const' (2/3) -l)/b) 
EH2=2*pi-EHl 

% Berechnung der Wellenzahl : 
integrand=-oniega. *w ; 
deltax=(xende-xanf ang) / anzahl; 
for j = 1:2 

integraK j , : )=cumsum(integrand( j , : ) ) *deltax; % Nummerische Integration 

end 

% Berechnung der ebenen Wellen: 

phased, :)= f aktor (1 ,:). * cos (integral (1 ,:)) ; 

phase(2,:)= f aktor (2 , : ) . * sin(integral(2, : ) ) ; 

neuphased , : )= faktor(l,:).* real(exp(i*omega*integral(l, : ))) ; 

% Berechnung der y-Werte fiir die Geraden durch den Ereignishorizont 
g=linspace (-2000 , 2000 , 1000) ; 
l=linspace(-5,5, 1000) ; 

% Plotten der Funktionen: 

7o Schall- und Hintergrundgeschwindigkeit 

subplot (3, 1,1) 

plot (x, sqrt (dichte_s) , 'b' ,x,v_s, 'r' ,EH1,1, 'k' ,EH2,1, 'k') 
axis ( [xanf ang xende 2]) 
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"/o Wellenzahl fiir beide Wellen 
subplot (3, 1,2) 

plot(x, integrand ( 1, :), 'r' ,x, integrand(2, : ) , 'b' ,EHl,g, 'k' ,EH2,g, 'k') 
axis( [xanfang xende -2000 2000]) 

°/o Ebene Wellen 
subplot (3, 1,3) y„ 

plot (x, phased, :) , 'r' ,x,phase(2, :) , 'b' ,EH1,1, 'k' ,EH2,1, 'k') 

axis ( [xanfang xende -5 5]) 

% Vergrofierung eines Ausschnitts aus dem vorhergehenden Plot 
figure 

plot (x, phased, :) , 'r' ,x,phase(2, :) , 'b' ,EH1,1, 'k' ,EH2,1, 'k') 
axis([EHl-0.01 EHl+0.01 -2 2]) 



D.2 Programm zur Simulation eines Wellenpakets im ringformigen 
Kondensat 



°/o Festlegung der Konstanten 



signia=0 . 1 ; 

k0=3; 

X0=3; 



°/o St andardabwei Chung 
°/o Max. des Wellenpakets im Impulsraum (t=0) 
% Max. des Wellenpakets im Ortsraum (t=0) 



Xmin=2 ; 
Xmax=6 ; 
NX=1000; 
Kmin=0 ; 
Kmax=3 ; 
NK=NX ; 
Tmin=0 ; 



% Anzahl der Intervalle im Ortsraum 



°/o Anzahl der Intervalle im Impulsraum 
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Tmax=-19; 
NT=50; 



% Anzahl der Intervalle im Zeitraum 



% Hauptteil 

7o Diskretisierung der Raume: 
n=l:NX; 

Xi=Xiiiin+(n-l) / (NX-1) * (Xmax-Xmin) ; 
DX=(Xmax-Xmin)/(NX-l) ; 
Kj=Kmin+(n-l)/(NK-l)*(Kmax-Kmin) ; 
DK=(Kmax-Kniin)/(NK-l) ; 
t=l:NT; 

Tn=Tmin+(t-l)/(NT-l)*(Tmax-Tniiii) ; 
DT=(Tmax-Tniin)/(NT-l) ; 

7o Berechnung der ebenen Wellen (siehe vorheriges Kapitel) 
const=l ; 

vs=const/(l+b*cos(XO) ~2) ; 
rhos=(l+b*cos(XO) ) ; 
A=sqrt (abs (vs-(-vs+sqrt (rhos) ) ) ) ; 
a=l/ (-vs+sqrt (rhos) ) ; 

dichte_s=(l+b*cos(XO)) ; 
v_s=(const . /dichte_s) ; 
EHl=acos ( (const' (2/3) -1) /b) 
EH2=2*pi-EHl 

AX=sqrt (abs (v_s-(-v_s+sqrt (dichte_s) ) ) ) ; 
aX=l . / (-v_s+sqrt (dichte_s) ) ; 
FX=cumsuiii(aX) *DX ; 

% Fourier-Analysen: 

% Fourier-Koef f izienten im flachen Raum (siehe Gl.(9.35)) 
Psi=l/sqrt(sigma)*exp(-(Xi-XO) . "2/(2*sigma) ) . *exp(i*kO*Xi) ; 
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EE=[]; 
MM= [] ; 
GG=[]; 
for nn=l:NK 

7o Fourier-Koef f izienten im "normalen" Flachen im Implusraum (siehe Gl.(9.37)) 
EE=[EE;exp(-i*Kj (nn)*Xi*a)] ; 

% zur Vorbereitung der Fourier-Analyse im gekrummten Raum 

7o in einer Zeile sind die Frequenzen konstand, in einer Spalte die Wellenzahl 
GG=[GG;exp(i*Kj*FX(nn)) ./AX(nn)] ; 

end 



°/o Berechnung der Fourier-Koef f izienten im fast flachen Raum (siehe Gl.(9.37)) 
FTPsi=(exp(i*Kj*Xmin*a) . ') . * (A*a/pi*EE*Psi . 0*DX; 



M=moviein(NT) ; % Graf ische Darstellung 
for ct=l:NT 

% Fourier-Analyse in der Zeit (siehe Gl.(9.34)) 
Psixt=GG*(FTPsi.*exp(-i*Kj*Tn(ct)) . ') ; 
plot (Xi , abs (Psixt) . "2) 
M(ct)=getframe; 

end 
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Anhang E 

Losung der 

Bogoliubov-Gleichungen fiir das 
ringformige Kondensat durch 
Fourier- Transformation 



E.l Berechnung von /, 



np 



fnp = ^ /o d7e-^("-P)^c(7)2 Die Schallgeschwindigkeit ist (siehe Gl. lFmil l 



i- j^J d7e-'("-^)^c2(7) = ^^£^ d7e-^(""^')^ (1 + 6005(7)) 



±— / d7e-*("-P)T ± / (i7e-^("-P)^cos(7) 



27r 
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np 



E.2 Berechnung von 



2n 



Jq^ dje ^''^pvo{j) Die Hintergrundgeschwindigkeit ist Gl. (|9.14|) 



u(7j = m 
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(1 + 6cos(7) 



1 + 6 COS (7) 
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